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Section 0: Introduction 

This paper introduces a rigorous computer-assisted procedure for analyzing hyperbolic 
3-manifolds. This technique is used to complete the proof of several long-standing rigidity 
conjectures in 3-manifold theory as well as to provide a new lower bound for the volume 
of a closed orientable hyperbolic 3-manifold. 

Theorem 0.1: Let N be a closed hyperbolic 3-manifold. Then 

1) If /: M —>■ N is a homotopy equivalence where M is a closed irreducible 3-manifold, 
then / is homotopic to a homeomorphism. 

2) If /, g: M — *■ N are homotopic homeomorphisms, then / is isotopic to g. 

3) The space of hyperbolic metrics on iV is path connected. 

Remarks: Under the hypothesis that M is hyperbolic, conclusion i) follows from Mostow's 
rigidity theorem [Mo]. Under the hypothesis that N is Haken, conclusions i)-ii) follow from 
Waldhausen [W]. If iV is both Haken and hyperbolic, then iii) follows by combining [Mo] 
and [W] . Since non-Haken manifolds are necessarily orientable we will from now on assume 
that all manifolds under discussion are orientable. 

Theorem 0.1 with the added hypothesis that a closed geodesic 6 C N has a non 
coalescable insulator family was proven by Gabai (see [G]). Thus Theorem 0.1 follows from 
[G] and the main technical result of this paper which is, 

Theorem 0.2: If S is a shortest geodesic in a closed orientable hyperbolic 3-manifold, 
then S has a non-coalescable insulator family. 

Remarks: If 5 is the core of an embedded hyperbolic tube of radius (ln3)/2 = 0.549306... 
then S has a non-coalescable insulator family by Lemma 5.9 of [G]. In this paper we 
establish a second condidition, sufficient to guarantee the existence of a non-coalescable 
insulator family for 5. That is if Corona (S) < 2n/3. (Corona(5) < 27r/3 if tube radius(5) > 
(ln3)/2). We use the expression U N satisfies the insulator condition" when there is a 
geodesic 5 which has a non-coalescable insulator family. 

We prove Theorem 0.2 by first showing that all closed hyperbolic 3-manifolds, with 
seven families of exceptional cases, have embedded hyperbolic tubes of radius (ln3)/2 about 
their shortest geodesies. Conjecturally, up to isometry, there are exactly six exceptional 
manifolds (see Conjecture 4.2). Second, we show that any shortest geodesic 5 in six of the 
seven families has Corona(5) < 27r/3. Finally, we show that the seventh family corresponds 
to Vol3, the third smallest known hyperbolic 3-manifold, and that the insulator condition 
holds for Vol3. Each of the three parts of the proof is carried out with the assistance of a 
rigorous computer program. 

Here is a hint why Theorem 0.2 might be amenable to computer proof. If a shortest 
geodesic 5 in a hyperbolic 3-manifold N does not have a (ln3)/2 tube then there is a 
2-generator subgroup G of tti(N) = V which also does not have that property. That is, 
after identifying N = H 3 /r and letting Z = H 3 /67, then a shortest geodesic in Z does 
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not have a (ln3)/2 tube. G is a group of 2 generators, generated by / and w, where / e V 
is a primitive hyperbolic isometry whose fixed axis So C H 3 projects to 5, and w is an 
isometry of H 3 which takes do to a nearest translate. Here So is a lift of 5 to H 3 . 

The space of relevant 2-generator groups in Isom(H 3 ) naturally lives in C 3 . We show 
that except for seven small regions in C 3 the shortest geodesic in any discrete, torsion free, 
parabolic free 2-generator group must have a (ln3)/2 tube. Further, if Corona(5) > 27r/3, 
then there is a 2-generator subgroup with that property. That is, there is a 2-generator 
subgroup G of 7ii(N) such that if N\ = H 3 /G, then Corona(5i) = Corona(5) > 27r/3 for 
some shortest geodesic 5i in N±. We show that away from a single small open set in C 3 , 
every discrete, torsion free, parabolic free 2-generator group G satisfies Corona(5) < 27r/3, 
where 5 is a shortest geodesic in H 3 /C We finally show that the exceptional open set 
contains a unique manifold which is Vol3. In fact we show that if a shortest geodesic 5 in 
N satisfies Corona(5) > 27r/3, then N = Vol3. A variant of the above arguments shows 
that Vol3 satisfies the insulator condition. 

This paper is organized as follows. In Chapter 1 we describe a space V C C 3 which 
naturally parametrizes all relevant 2-generator groups. We explain how a theorem of 
Meyerhoff as well as elementary hyperbolic geometry considerations imply that we need 
only consider a compact portion V of C 3 . We explain in detail the plan for proving 
Theorem 0.2. We will actually be working in the parameter space W = exp('P). The 
technical reasons for working in W is described at the end of this section. In Chapter 2 we 
describe and prove the necessary results about the Corona function. In Chapter 3 we prove 
that the exceptional open set in C 3 contains only Vol3. Also if 8 is a shortest geodesic 
in a closed orientable hyperbolic 3-manifold N and Corona(5) > 2/t/3, then N =Vol3. 
Nonetheless, we are able to show that Vol3 satisfies the insulator condition. In Chapter 4, 
we prove some applications, one of which is discussed briefly below. 

In Chapters 5 through 8 we address the computer-related aspects of the proof. Here, 
the method for describing the decomposition of the parameter space W into sub-regions is 
given, and the conditions used to eliminate all but seven of the sub-regions are discussed. 
At the end of this chapter, the first part of a detailed example is given. Eliminating a 
sub-region requires that a certain function is shown to be bounded appropriately over the 
entire sub-region. This is carried out by using a first-order Taylor series approximation 
of the function together with a remainder bound. Our computer version of such Taylor 
series with remainder bounds is called an AffApprox and in Chapter 6, the relevant theory 
is developed. At this point, the detailed example of Chapter 5 can be completed. 

Finally, in Chapters 7 and 8, round-off error analysis appropriate to our set-up is 
introduced. Specifically, in Chapter 8, round-off error is incorporated into the AffApprox 
formulas introduced in Chapter 6. The proofs here require an analysis of round-off error 
for complex numbers, which is carried out in Chapter 7. 

We used two rigorous computer programs in our proofs — verify and fudging. These 
programs are provided in Appendices 1 and 2. Actually, fudging is a variation of verify and 
as such we only provide the sections of fudging that are changed. The proofs amount to 
having verify and fudging analyze several computer files. These computer files are available 
from the Geometry Center. Details about how to get them and the programs can be found 
at 
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http : //www . geom . umn . edu : /locate/HomotopyHyperbolic 

A consequence of this work is that either the shortest geodesic in a closed orientable 
3-manifold N has a 1.059191579962 . . . /2 tube or N = V613. The volume of Vol3 is 1.01... 
and by [GM2] if N has a log(3)/2 = 0.549306 . . . tube about a geodesic, then the volume 

of N is greater than 0.16668 This leaves some exceptional cases that can be analyzed 

using data provided by verify. In any case, we obtain 

Theorem 4.5: If A is a closed hyperbolic 3-manifold, then vol(A) > 0.16668 

Remark: The best published lower bound for volume is 0.001 by [GM1], which improved 
the lower bound of 0.0008 of [M2]. 
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Chapter 1: Killer Words and the Parameter Space 

Definition 1.1: We will work in the upper-half-space model for hyperbolic 3-space. All 
isometries will be orientation preserving. If / is an isometry, then we define Relength(/) = 
inf{p(x, f(x)) | x G H 3 }. Thus Relength(/) = if and only if / is either a parabolic or 
elliptic isometry. If Relength(/) ^ 0, then / is hyperbolic and fixes a unique geodesic a in 
H 3 . In that case a is oriented (the negative end being the repelling fixed point on S^) and 
the isometry / is the composition of a rotation of t (mod 2tc) radians along a (the sign of 
the angle of rotation is determined by the right-hand rule) followed by a pure translation 
of H 3 along a of / = Relength(/). We define length(/) = l + it. If a is an oriented geodesic 
in H 3 , then it makes sense to talk about an I + it translation of H 3 along a, even when 
/ < 0. 

If / is elliptic, then / is a rotation of t radians about some oriented geodesic. If / is 
elliptic, we define length(/) = \t\i, the absolute value accounting for the arbitrariness of 
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the orientation of the fixed geodesic. If / is parabolic or the identity, we define length(/) = 
+ iO. So, for all isometries we have that Relength = Re(length). 

Definition 1.2: If G is a subgroup of Isom(H 3 ), then we say that / is an element of 
smallest length in G if Relength(/) < Relength((7) for all g E G, g ^ id. 

Convention 1.3: Let B denote the oriented geodesic £(0, 0, 1), with negative end (0, 0, 0). 
Let C denote the oriented geodesic with negative endpoint (—1, 0, 0) and positive endpoint 



Lemma 1.4: If the isometry / is represented by the matrix A E PSL(2, C), then 



where the branch of Arccosh with positive real values is taken, unless the real part is zero 
in which case the non-negative imaginary part is taken. 

Proof: Because trace is a conjugacy invariant, we can normalize our set-up via conjugation 
and assume that the axis of / is B. As such, A is a diagonal matrix, with diagonal entries 
p and p~ x . 

The action of A on the bounding complex plane is simply multiplication by p 2 . Extend- 
ing this action to upper-half-space in the natural way rotates the z— axis by angle arg(p 2 ) 
and sends (0,0,1) to (0, 0, \p\ 2 )- Thus, Im(length(/)) = arg(p 2 ) = Im(ln(p 2 )) and, using 
the hyperbolic metric, Re(length(/)) = ln(|p| 2 ) = Re(ln(p 2 )). That is, length(/) = ln(p 2 ). 
Thus, we need only show that 2 ln(p) = 2Arccosh(trace(A)/2). But this follows because 
cosh(ln(p)) = (p + p~ 1 )/2 = trace(A)/2. | 

Definition 1.5: If a, r are disjoint oriented geodesies in H 3 which do not meet at 
infinity, then define distance(a, r) = length(w) where w E Isom(H 3 ) is the hyperbolic 
element which translates H 3 along the unique common perpendicular between a and r 
and which takes the oriented geodesic a to the oriented geodesic r. The oriented common 
perpendicular from a to r is called the orthocurve between a and r. The ortholine between 
a and r is the complete oriented geodesic in H 3 which contains the orthocurve between a 
and t. 

If a, t intersect at one point in H 3 then slight changes must be made in the above 
definition. The ortholine has no natural orientation, the orthocurve is a point, and w is 
an elliptic isometry. 

If a, t intersect at infinity, then there is no unique common perpendicular, hence no 
ortholine, and distance(cr, r) = + z0, or O + in depending on whether or not a and r point 
in the same direction at their intersection point (s) at infinity. 

Lemma 1.6: distance(cr, r) = distance (r, a) I 

Lemma 1.7: If the isometry / is represented by the matrix 



(1,0,0). 



length(/) = 2Arccosh(trace(A)/2), 



A = 




E PSL(2,C), 
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then 2distance(/(S), B) = 2Arccosh(a<i + be). Again, the branch of Arccosh with positive 
real values is taken, unless the real part is zero in which case the non-negative imaginary 
part is taken. 

Proof: In the case where B and f(B) do not intersect at infinity, we will compute the 
length of h, the square of the transformation taking B to f(B) along their ortholine. 
h = (/ ° P ° ° P where p is 180-degree rotation about B and hence (/ o p o is 
180-degree rotation about f(B). p and / are represented by the matrices 

Hence, h = (/ o p o o p can be computed to have representation 

/ ad + bc 2ab \ 
\ led ad + be J 

Using Lemma 1.3, we have that 

2distance(/(S), B) = length(/i) = 2Arccosh(trace(/i)/2) = 2Arccosh(a<i + be). 

If / fixes the point (0, 0, 0) at infinity, then c = 0, ad = 1 and the formula holds. 
Similarly for the other cases in which B and f(B) intersect at infinity. 
Note that this formula only defines distance(/(S), B) modulo ztt, but this is sufficient for 
our purposes. I 

Definition 1.8: Let 5 be a geodesic in the hyperbolic 3-manifold N. Then tuberadius(d) = 
sup{r | there exists in N an embedded D 2 x S 1 of radius r centered about the geodesic 8}. 

Lemma 1.9: Let 5 be a geodesic in the hyperbolic 3-manifold N and {5i}i>o be the set 
of its distinct lifts to H 3 , then tuberadius(5) = min{Redistance(5o 5 Si) I i 7^ 0}. I 

Definition 1.10: We define an open subset V of C 3 which naturally parametrizes the 
collection of conjugacy classes of 2-generator subgroups G C Isom(H 3 ) with specified 
generators /, w where / is hyperbolic and w is not parabolic. By conjugating G we can 
assume that / is a positive translation of H 3 along the geodesic B, and that the orthocurve 
from w~ 1 (B) to B lies on C on the negative side of C n B. Associated to {67, /, w}, the 
group G with specified generators /, to, is the parameter (L, D, R) = (/ + it, d + ib,r + ia) 
where I + it = length(/) and d + ib = distance(w(S), B). The complex number r + ia 
is defined as follows. The isometry w is the composition of two isometries, first a d + ib 
translation of H 3 along C, which takes w~ 1 (B) to £?, as oriented geodesies, followed by 
an r + ia translation along B. 

We are primarily interested in the set T C C 3 which parametrizes all conjugacy 
classes of triples {67, /, w} where G is a group generated by a shortest element / which 
(positively) translates B and w G G takes B to a nearest translate w(B) such that 
— Relength(/)/2 < Redistance( (ortholine from w~ 1 (B) to £?), (ortholine from B to 
w(B))) < Relength(/)/2. By conjugating G we can assume both that / is a positive 
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Remark 1.11: Said another way, T corresponds to those parameters such that / is the 
real length of a shortest element of G, d is the real distance between B and a nearest 
translate, and — 1/2 < r < 1/2. In what follows, it is essential to remember that an element 
a of V corresponds not only to a group G, but a group with two special generators. When 
a G T', then two (i.e. / and d) of a's six real parameters correspond to invariants of 
{G,f,w}. 

We are only interested in the subset of T corresponding to parameters a with d < 
ln(3). The following two propositions imply this subset of T lives in a compact subset of 
V. 

Proposition 1.12: All closed geodesies of length less than 0.0979 in all hyperbolic 3- 
manifolds have embedded solid tube neighborhoods of radius In 3/2. 

Proof: In [Ml] it is proven that a geodesic of length x + iy has an embedded solid-tube 
neighborhood of radius r(x + iy) satisfying 



sinh 2 (r(;r + iy)) = max -(— — —. ^ 1) where k(x,y,n) = cosh(ra) — cos(ra/). 

nez + 2 k(x, y,n) 

Of course, we restrict to x + iy values which produce positive radii r(x + iy) by means of 
this formula. It is easy to compute that for a given x + iy we need to have n for which 
< k(x, y, n) < —1 + \/2 to produce a positive radius tube by this method. 

The function \{ ^\ 2k — 1) is decreasing on the interval (0, — 1 + y/2). It is therefore 
easy to solve for the range of k values that produce radii r greater than In 3/2. In fact, 
positive k less than 0.3397 work. 
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Thus, to complete the proof of this proposition, we need to show that when a geodesic 
has real length x less than 0.0979, that for all angles y there exists a positive integer n 
for which k(x,y,n) is less than 0.3397. Because cosh is an increasing function, we can 
restrict our analysis to x = 0.0979. Thus, we need only show that given any angle y, we 
can find a positive integer n such that cosh(n0.0979) — cos(ny) < 0.3397. When n > 8 we 
can compute that cosh(n0.0979) — cos(ny) > 0.3397, and we therefore restrict to positive 
integers n < 8. 

We now consider angles y. Because cos is an even function, we need only consider 
y G [0, 7r]. We will cover [0, n] by 11 overlapping closed sub-intervals Oi each of which has 
an associated positive integer rii for which cosh(ni0.0979) — cos(riiy) < 0.3397 is true for 
all y e Uj. 



cr = 


[0.000, 0.843] 


n 


= 1 


(7l = 


[0.835,0.960] 


m 


= 7 


02 = 


[0.951,1.143] 


n 2 


= 6 


3 = 


[1.123,1.391] 


n 3 


= 5 


(T4 = 


[1.386, 1.755] 


714 


= 4 


(T5 = 


[1.733, 1.858] 


n 5 


= 7 


cr 6 = 


[1.832,2.357] 


n 6 


= 3 


0-7 = 


[2.334,2.3792] 


n 7 


= 8 


(7 8 = 


[2.3789,2.647] 


n 8 


= 5 


cr 9 = 


[2.630, 2.755] 


n 9 


= 7 


O"io 


= [2.730, ?r] n 


10 = 


= 2 



Proposition 1.13: If the shortest geodesic in a closed hyperbolic 3-manifold has length 
greater than or equal to 1.29, then it has an embedded solid tube neighborhood of radius 
In 3/2. 

Proof: Consider the following folklore result: The shortest geodesic in a closed hyperbolic 
3-manifold has an embedded solid tube neighborhood of radius 1/4 where / is the real 
length of the shortest geodesic. The proof is simple: Expand a solid tube around the 
shortest geodesic, if it hits itself before a radius of 1/4 then we will construct a loop of 
length less than /, a contradiction to "shortestness." Drop the two obvious perpendiculars 
from the hitting point down to the core geodesic. Consider the following loop — down one 
perpendicular, follow the shorter direction on the core geodesic, up the other perpendicular. 
This non-trivial loop has length less than 1/4 + 1/2 + 1/4 < I. 

We now improve on this loop. Replace the first half of the journey by the hypotenuse 
of the right triangle formed by the first perpendicular and the first half of the shorter arc 
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along the core geodesic. Replace the second half of the journey by the hypotenuse of the 
right triangle formed by the second perpendicular and the second half of the shorter arc 
along the core geodesic. Using the hyperbolic Pythagorean Theorem (see [F]) coshc = 
(cosh o) (cosh b) with a = In 3/2 and 6 = 1/4 we get that the length of the constructed loop 
is 2 cosh - 1 (cosh(ln 3/2) (cosh(//4))) and this is less than / when / > 1.29, by a calculation 
and the fact that 2 cosh _1 (cosh(ln3/2)(cosh(//4)) — / is a decreasing function of /. 
In fact, we can solve explicitly for the value of / at which 

2cosh" 1 (cosh(ln3/2)(cosh(//4)) - / = 0. 

Noting that cosh(ln3/2) = we get -^(cosh(//4))) = cosh(//2). Using a half-angle 
formula for cosh(//2) we get ^(cosh(//4))) = 2cosh 2 (//4) - 1. Setting x = cosh 2 (//4) we 
get the quadratic -^=x = 2x 2 — 1. Solving and substituting, we get / = 1.289784... | 

Definition 1.14: Let V C V be those parameters a = (I + it,d + ib,r + id) such that 

a) -7T < t < 

b) < r < 1/2 

c) 0.0978 < I < 1.29 

d) 1/2 <d< ln(3) 
Define T = T'nV. 

Lemma 1.15: If a = (I + it, d + ib, r + ia) G T' has d < ln(3) and corresponds to a 2- 
generator group {G a , f a , w a }, then there exists a parameter (3 G T with associated group 
{G/3, fp,wp} such that Gp is conjugate to G a . 

Proof: d < 1/2 is eliminated from consideration by the first paragraph of the proof of 
Proposition 1.13, and the definition of T . If d < Zn(3), then 0.0978 < / < 1.29 by Propo- 
sition 1.12. If for the triple {G, /, w} we have — 1/2 < r < 0, then the triple {G, /, is 
conjugate to an element of T whose new r-parameter is — r. Thus we can assume that b), 
c), and d) hold for the relevant {G, f, w}. 

This leaves property a). Conjugating G by a reflection in the geodesic plane spanned 
by B and C changes the t-parameter to —t (mod 2ir). The effect on b and a is irrelevant. 
I 

By [G; Lemma 5.9] a closed orientable hyperbolic 3-manifold iV satisfies the insulator 
condition provided that tuberadius(5) > ln(3)/2 for some geodesic 5 C N. Thus we are led 
to ask 

Question 1.16: List all closed orientable hyperbolic 3-manifolds iV possessing a shortest 
geodesic 5 such that tuberadius(5) < ln(3)/2. 

Remarks 1.17: i) Using a J. Weeks-modified version of Snappea it was known exper- 
imentally that any shortest geodesic in Vol3 has a .415... tube (see [G]). Conjecturally, 
up to isometry, there are a total of six manifolds in the list answering Question 1.16 (see 
Theorem 1.30 and Remark 1.31 iii, and Theorem 4.xx). 

ii) If a shortest geodesic 5 in N satisfies tuberadius(5) < ln(3)/2, then iV can be 
expressed as H 3 /r where a lift of 5 is the geodesic B, and C is an ortholine between B 
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and a nearest T— translate. Thus N gives rise to an element a G T. In fact N may give 
rise to finitely many different elements of T. Thus we need to investigate. 

Question 1.18: Name all parameters a = (I + it, d + ib, r + ia) G T. 

We now describe our method of (partially) answering Question 1.18. There is a 
technical point to mention: starting with Definition 1.22, we will work in the space W D 
exp('P), but for now we will describe the results in terms of the unexponentiated space V . 

We will partition V into about one billion regions {Vi} and show that T is disjoint 
from all but seven small such regions. Suppose that Vi is a region of this partition and 
a G Vi. Let h be a word in the letters /, w and their inverses. Associated to the parameter 
« = (la + ita, d a + ib a , ra + iaa) there are the group elements f a , w a and hence h a . Suppose 
that h a ^ f™. We ask 

a) Is Relength(/i a ) < Relength(/ a ) = / a ? 

b) Is Redistance(/i a (S), B) <Redistance(w a (S), B) = d a ? 

If either a) or b) is true, then a ^ T. 

Now let (3 G Vi, with fp,wp, and hp the associated hyperbolic isometries. If say a) 
is true for a then so is the statement Relength(/ig) < Relength(/ ( g) = lp for (3 sufficiently 
close to a. Thus we can show that T D Vi = if we can find an a for which say statement 
a) is true, and then use first-order Taylor approximation (with error /remainder term) to 
show that the corresponding statement holds for all f3 G Vi. 

Definition 1.19: A word h in w, /, ty -1 , f~ l for which statement a) (resp. b)) holds 
(non-trivially) for each (3 G Vi is called a killer word for V% with respect to contradiction 
a) (resp. b)). 

Summary 1.20: With seven exceptions, to each of the approximately one billion regions 
partitioning V, we will associate a killerword and a contradiction. 

Remark 1.21: Computers are well suited for partitioning a region such as V into many 
sub- regions {Vi}, and finding a killerword hi which eliminates aij G Vi due to contradiction 
Cj. Depending on the contradiction, we find computable expressions for approximations of 
the values of Relengtt^/i^) or Redistance(/i ( g(i?), B) and thus use the computer to eliminate 
all of Vi . 

There are a number of difficulties in executing this procedure. First, a uniform mesh 
of the partition would yield far too many sub-regions to be handled by computer. In fact 
with 6 real parameters, refining a given mesh by a factor of 10 would change the partition 
size by a factor of 10 6 . Our method for refining the parameter space and the way the 
computer keeps track of the refinements are discussed briefly in Remark 1.27 and in more 
detail in Chapter 5. 

A second difficulty is finding the killerwords. In practice, most of the parameter space 
is eliminated by killerwords of length less than 7, but a number of spots need killerwords of 
length 10 and a few regions need killerwords of length 35. A brute force enumeration and 
testing of the various words would take far too long. Note that there are more than 70,000 
words of length 10 and 4 x 3 29 words of length 30. Techniques for finding killer words can 
be found in [Txxx]. 
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Finally, there is the issue of rigor. The main difficulty in making the plan work 
rigorously is that we need to bound the difference between what the computer thinks a 
value is and what the value actually is. In particular we need to control roundoff error, 
which becomes quite significant when one does a large number of multiplications, e.g. 
in the computation of length(/i a ) when h is a 35 letter word. Another issue is to make 
sure that hp is bounded away from f™, and in particular bounded away from id, when 
computing say Re\ength(hp). A large portion of this paper is devoted to addressing these 
issues. See Remark 1.30 for a more detailed discussion. 

Definition 1.22: Let 

W = {(x ,x 1 ,x 2 ,x 3 ,x 4 ,x 5 ) : \xi\ < 4 x 2 (5_i)/6 for i = 0,1,2,3,4,5} 

D exp(P) = {(xo,Xi,X2,X3,X4,X 5 ) I 

x + ix 3 = exp(e), x\ + zx 4 = exp(/), x 2 + ix 5 = exp((7)where(e, /, g) E V} 



L' = exp(L) = exp(7 + it), D' = exp(D) = exp(d + ib), R' = exp(R) = exp(r + id). 

Remarks 1.23: i) We work with W instead of exp('P) because we want our initial region 
to be a (6-dimensional) box that is easily sub-divided. This has the side-effect that certain 
sub-boxes Wj of W will be eliminated because they are outside of exp('P) rather than by 
the analogues of conditions a) and b) above. The entire collection of conditions is given in 
Section 5. 

ii) All the ideas expressed in 1.18-1.21 will be carried out in the parameter space W 
rather than the space V. That is mainly because of Lemmas 1.24 - 1.26 which demon- 
strate that while working in W one need only understand the basic arithmetic operations 



hi) The reason for choosing the co-ordinates of W so that V = xq + ixs, D' = 
x\ + ix^, R' = X2 + ix^ was to gain a mild computer advantage. 

Lemma 1.24: If (L', D' ', R') e W and /, w are the generators of the associated group G 
then 



and let 



S = exp(T). 



Also, let 



+,-, x,/ 5v /. 



a) 




b) 




where ch = ( 




+ V 




) /2 and sh = ( 




V 
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Proof: a) By our set-up we have that the (oriented) axis of / is B. Following the proof 
of Lemma 1.3, 

Matrix[/] =(g ° /p ) 

where p = exp(length(/)/2) = A/exp(length(/)) = ^exp(L) = \fU 

b) w = (3 o a where j3 is translation of distance R along B and a is translation of 
distance D along C. Thus, 

Matrix^] =(^f 

and Matrix [a] can be computed to be 

/cosh(D/2) sinh( J D/2)\ 
^ sinh(D/2) cosh(D/2) J ' 

But cosh( J D/2) = (exp(D/2) + exp(-D/2))/2 = (Vd 7 + l/VD')/2 = eh and similarly for 
sh. Thus, 

Matrix[a]=(^ ^) 

and b) follows by matrix multiplication. | 

Lemma 1.25: If h e Isom(H 3 ) is represented by the matrix 

A ={l d)^SL(2,C), 

then 

a) exp(Relength(/i)) = |trace(A) ± v /(trace(A)/2) 2 - 1| 2 

b) exp(Redistance(/i(B), B)) = |orthotrace(A) ± A/(orthotrace(A)/2) 2 - 1| 

where orthotrace(A) = ad + be. 
Here the +, — produce reciprocal values for exp(Relength(/i)), and we take the one pro- 
ducing the larger value, unless the value is 1, in which case there is no choice. 

Proof: Because cosh(x) = (exp(x) + exp(— x))/2 it is easy to compute that cosh _1 (x) = 
log(x±y/x 2 - 1). Of course, x — \Jx 2 — 1 and x + \Jx 2 — 1 are inverses, which corresponds 
to the fact that cosh _1 (x) for x ^ 1 consists of values differing by a factor of —1. 

a) exp(Relength(/i)) = | exp(length(/i))| = | exp(2Arccosh(trace(A)/2))| = |(trace(A)/2) ± 
a/ (trace(A)/2) 2 — 1| 2 where the second equality follows from Lemma 1.2. 

b) exp(Redistance(/i(S), B)) = | exp(distance(/i(S), B))\ = | exp(Arccosh(a<i + bc))\ = 
|(aei-|- be) ± a/ (ad + be) 2 — 1| where the second equality follows from Lemma 1.7. I 

Remarks 1.26: i) It follows from Lemma 1.25 that if h is a word in /, w and their inverses, 
then for any parameter value a G W, exp(Relength(/i a )), and exp(Redistance(/i a (S), B)) 
can be computed using only the operations +, — , x, /, ^j. | 
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ii) During the course of the computer work needed to prove the main theorems, the 
parameter space W was decomposed into sub-boxes by computer via a recursive sub- 
division process: Given a sub-box that is being analyzed, either it can be handled, or 
it cannot. If it cannot be handled, it is sub-divided in half by a hyper-plane {xi = c} 
(where i runs through the various co-ordinate dimensions cyclically) and the two pieces 
are analyzed separately. And so on. 

As such, a sub-box of W can be described by a sequence of O's and l's where means 
"take the lesser Xi values" and 1 means "take the greater xi values." Remarkably, for the 
decomposition of W into sub-boxes, all the sub- box descriptions could be neatly encoded 
into one tree (although in practice we found it preferable to use several trees to describe 
the entire decomposition). This is described in Chapter 5. 

iii) In the following proposition, seven exceptional boxes are described as sequences 
of O's and l's. Four of the boxes — Xq, X4, X5, A^ — are each the union of two abutting 
sub-boxes, X = X 0a U A & and so on. It is a pleasant exercise to work through the fact 
that they abut. 

It is also a pleasant exercise to calculate by hand the co-ordinate ranges of the various 
sub-boxes. For example, the range of the last co-ordinate (i.e., x 5 ) of the sub-box 

x 6a = 111000000001000111 111111110101001111 011111010111111111 

110001001011000111 

is found by taking the 6th entry, the 12th, entry, the 18th entry and so on. These entries 
are 011111111111. The first entry (0) means take the lesser x 5 values, and produces the 
interval [—4,0]. The second entry (1) means take the greater x$ values, and produces the 
interval [—2,0]. The third entry (1) produces [—1,0]. Continuing, we see that X^ a has 
— 2 -9 < £5 = Im(R') < 0. The other co-ordinates can be computed in the same fashion, 
although they must at the end be multiplied by the facter 2^ 5- ^/ 6 (see the definition of 
the initial box W). The range of co-ordinate values is given for each of the seven boxes. 
Finally, two quasi-relators are given for each sub-box A , Xi, . . . , X G . 

Definition 1.27: A quasi-relator in a sub-box X of W is a word in /, that 
is close to the identity throughout X and experimentally appears to be converging to the 
identity at some point in X. In particular, a quasi-relator rigorously has Relength less than 
that of / at all points in X. 

Proposition 1.28: SD(W— U n =i 1 ^-n) = where the X n are the exceptional sub-boxes 

X Qa = 001000110111110001 101001010101011001 011011010111101101 

100001101101000111 010001110101100101 1101110111110100 

X Qb = 001001110110110000 101000010100011000 011010010110101100 

100000101100000110 010000110100100100 1101100111100100 
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Xn = 



( -0.840655162503.. 
-0.840642408899.. 

0.999979499517.. 
-2.137267196028.. 
-2.137295441962 . . 
V -0.000061035156. 



. < Re(L') < -0.840600786360 . 
. < Re(D') < -0.840593965263. 
. < Re(R') < 1.000022657890. . 
. < Im(L') < -2.137228746289. 
. < Im(D') < -2.137226932315 
. . < Im(R') < 0.000061035156 . 



X quasi-relators: 

r\ = fwFwwFwfww 
T2 = FwfwfWfwfw 



X x = 001000110001110110 011101000110111110 100010110000100011 

101101001101001000 110101011000000100 000 



Xi = 



/ -1.348528333122.. 
-0.543343817104 . . 

0.903908961497.. 
-2.661029541660. . 
-2.858770529287.. 
V -1.471679687500.. 



. < Re(L') < -1.348310828552 . 
. < Re(D') < -0.543150042561 . 
. < Re(R') < 0.904081594988 . . 
. < Im(L') < -2.660721943747. 
. < Im(D') < -2.858496490701 
. < Im(R') < -1.471435546875 . 



X\ quasi-relators: 

n = FFwFWFWfWFWFwFFww 
T2 = FFwwFwfwfWfwfwFww 



X 2 = 001000110101010010 101010110001100101 110111100001101010 

111100100000010001 111100 



Xo 



( -1.787017545957 
-1.074286063490. 
0.741633479486. 
' -2.272533378081. 
-2.718462773249 . 
V -1.529296875000. 

quasi-relators: 

r\ = FwfwfWffWfwfwFww 
T2 = FFwFFwwFwfwfwFww 



. < Re(L') < -1.785277509398 . . 

< Re(D') < -1.072735867150. . 

< Re(R') < 0.743014547418 . . . 

< Im(L') < -2.271302986431 . . 

< Im(D') < -2.717366618905. . 

< Im(R') < -1.528320312500 . . 



\ 



X 3 = 111000000001000110 011011101101011000 111101011110001100 

111111100110110000 0000100010100010 
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X* = 



( 0.581172210661 . . . < Re{V) < 0.581607219801 . . . 
1.156446469500 . . . < Re(D') < 1.156834018585 . . . 
1.404200819866 . . . < Re(R') < 1.404546086849 . . . 
-3.312214322575 . . . < Im(L') < -3.311906724662 . . 
-2.756280098119 . . . < Im(D') < -2.755732020947 . 
V -1.179687500000 . . . < Im(R') < -1.179199218750 . . 

quasi-relators: 

n = FFwfwFFwwFWFwFWfWFWffWFWfWFwFWFww 
r 2 = FFwfwFwfWfwfWWfwfWfwFwfwFFwwFWFww 



\ 



X 4a = 111000000001000110 011001001111101010 011110110110111101 

100011111110110110 10000111101 



Xi b = 111000000001000110 011001001111101010 111110010110011101 

000011011110010110 00000101101 



Xa 



x 4 



( 0.333217001023 . . . < Re(L' 
0.977398792251 . . . < Re(D' 
1.354137107330... 
-3.319596672476 . . . 
-2.825337821794 
V -1.225585937500 

quasi-relators: 

n = FFwfwFwfWfwfWfwFwfwFFwwFWFwFWFww 
r 2 = FFwfwFwfwFFwwFWFwFWfWFWfWFwFWFww 



< 0.334957037582 . . . 

< 0.978173890421 . . 

< 1.354827641296. . 

< -3.318981476651 . 

< -2.824789744622 



< Re(R') 

< Im(L') 

< Im{D') 

< Im(R') < -1.224609375000 . . . / 



X 5a = 001000110001110111 001111000101111111 101111100111001111 

000001111011110111 1 



X 5b = 001001110000110110 001110000100111110 101110100110001110 

000000111010110110 1 



X* = 



( -1.379848991182.. 

-1.379674742433.. 
0.999893182771 . . 

-2.537067582893 . . 

-2.536501152136.. 
V -0.001953125000. 



. < Re(L') < -1.378108954623. 
. < Re(D') < -1.376574349753. 
. < Re(R') < 1.002655318635 . . 
. < Im(L') < -2.534606799593 . 
. < Im(D') < -2.534308843448 
. . < Im(R') < 0.001953125000 . 
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X 5 quasi-relators: 

r\ = FwFWFwFwfwfWfwfw 
r 2 = FwfwfWfWFWfWfwfw 



X 6a = 111000000001000111 111111110101001111 011111010111111111 

110001001011000111 



X 6b = 111001000000000110 111110110100001110 011110010110111110 

110000001010000110 



( 1.378108954623 . . . < Re(L') < 1.379848991182 ... \ 
1.376574349753 . . . < Re(D') < 1.379674742433 ... 
0.999893182771 . . . < Re(R') < 1.002655318635 . . . 
-2.537067582893 . . . < Im(L') < -2.534606799593 . . . 
-2.536501152136 . . . < Im(D') < -2.534308843448 . . 
V -0.001953125000 . . . < Im(R') < 0.001953125000 . . . / 



Xq quasi-relators: 

n = FWFwFWfWFwFWFwfw 
T2 = FWFwfwFwfWfwFwfw 

Proof: Two computer files contain the data needed for the proof. The first computer file 
describes the partition of W into sub-boxes, and the second describes the killerwords and 
contradictions associated with each sub-box (other than the Xi). We have a computer 
program verify which shows that the killerwords in question actually do kill off their 
associated sub- boxes. This computer program addresses the issues of Remark 1.21. The 
code for verify is given in Appendix 1, although we encourage readers to produce their 
own verification programs. 

In addition, verify showed that the listed words were quasi-relators for the given sub- 
boxes. I 

Corollary 1.29: If S is a shortest geodesic in N, a closed orientable hyperbolic 3-manifold, 
then 

i) either tuberadius(5) > ln(3)/2 or exp( length (5)) G C(Xk) for some k e 0, . . . , 6 
where C(Xk) denotes the range of V values in the sub-box X^. | 

ii) either tuberadius(5) > ln(3)/2 or tuberadius(5) = Re(L>)/2 where exp(D) e V(Xk) 
for some k e 0, . . . , 6 where V(Xk) denotes the range of D' values in the sub-box X\~. | 

Experimental Theorem 1.30: Associated to each of the sub-boxes Xq, Xi, . . . , Xq of W 
is a closed orientable 3-manifold of Heegaard genus 2 with fundamental group generated 
by f,w with 2 relators r\, r 2 as given in Proposition 1.28. 

Experimental Proof: Experimentally, at some point in the sub-box Xi under consid- 
eration the quasi-relators are actually relators. Applying Berge [Be], it follows that the 
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2-generator 2-relator presentation < /, w; r±, > is the group presentation of a Heegaard 
genus 2 closed orientable 3-manifold. I 

Remarks 1.31: i) In Chapter 3 it is proven rigorously that Vol3, the hyperbolic 3-manifold 
with the third smallest volume known, is the unique manifold associated with sub-box Xq. 
As of this writing its not known if the other boxes have unique manifolds. 

ii) Except for Vol3, Riley's program POINCARE was the first to show (experimentally) 
that there is a closed orientable (hyperbolic) 3-manifold associated to each box. It provided 
a group presentation, which presumably could have been shown to be those of the above. 

hi) Experimental evidence suggests that the manifolds associated to and Xq are 
isometric hyperbolic 3-manifolds, isometric to the Weeks census manifold s479(— 3, 1). Also 
experimentally, the manifold associated with X<i is s778(— 3, 1) and the manifold associated 
with X 1 is u2678(2,l). 

iv) Berge [Be] provides the explicit Heegaard genus 2 diagrams for each of the mani- 
folds in Experimental Theorem 1.30, so that the dilligent reader can recover Dehn surgery 
descriptions of these manifolds. 

Remark 1.32: It might be best to say that the list of killerwords and contradictions used 
to prove Proposition 1.28 was generated in an artistic rather than systematic mathematical 
way. Nevertheless, to have a rigorous mathematical proof one need only prove that these 
words work. 

Remark 1.33: Recall that we use first-order Taylor approximations (with Remainder 
term) which we denote AffApprox's to show that a killer word which eliminates a point 
x G Wi eliminates all of Wi. In setting up a Taylor approximation system we avoided a 
considerable amount of calculation by using the parameter space W instead of working 
with V. 

In the parameter space W, all functions analyzed via Taylor approximations are built 
up from the operations +, — , x , /, . We prove combination formulas for these functions, 
which show how the Taylor approximations (including the Remainder term) change when 
one of these operations is applied to two AffApprox's. This is carried out in Chapter 6. 

To ensure that all of our computer calculations are rigorous, we use a round-off error 
analysis. Typically, this is done by using interval arithmetic on floating-point numbers. 
This is too slow for our purposes, and so we introduce round-off error at the level of Af- 
fApprox's and incorporate the round-off error into the Remainder term. This also requires 
developing round-off error for complex numbers. This is all done in Chapters 7 and 8. 

Chapter 2: The Corona Insulator Family 

The upshot of Proposition 1.28 is that if a closed orientable hyperbolic 3-manifold 
has a shortest geodesic which does not have an embedded In 3/2 tube then the parame- 
ters for its associated 2-generator subgroup(s) (G, /, w) must be in one of the sub-boxes 
Xq,Xi,...,Xq listed in that proposition. Nonetheless — as we shall see in this Chapter 
and the next — such manifolds have non-coalescable insulator families about their shortest 
geodesies. However, they might not be Dirichlet insulator families. 
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In this Chapter, we describe a new insulator family {/%} called the Corona insulator, 
and we describe a condition sufficient for this family to be non-coalescable — a condition 
which is weaker than the tuberadius(5) > ln(3)/2 sufficient condition for the Dirichlet 
insulator family. 

The reason Dirichlet insulator families for geodesies with solid tubes of radius greater 
than In 3/2 are non-coalescable is that the amount of visual angle taken up by the various 
insulators is less than 120 degrees, and thus there is no chance for tri-linking to occur. 
The visual angle (measured at one axis) for a member of the Dirichlet insulator family 
associated to two axes depends only on the real distance between the two axes. 

In contrast, the visual angle for a member of the Corona insulator family associated 
to two axes depends on the (complex) distance between the two axes. We now give this 
function, C, and name it the Corona function. After that we give a precise definition of 
the visual angle function, and prove that the Corona function is the proper visual angle 
function for the Corona insulator family. 

Definition 2.1: Let C : (0, oo) x S 1 — > R be defined by 



C(u,v) = |Im(Arccosh(l 



1 ± cosh(tt + iv) 

where ± is positive for — tt/2 < v < tt/2 and negative otherwise. 



-t))I 




In the following definition, it is helgja4otai3B§$J#&rttk& 
in the upper- half-space model of H 3 . 

Figure 2.1 



1.05 1.10 

afuiSiJaeing the z-axis 
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Definition 2.2: If a C H 3 is a geodesic, then S^—da is parametrized by S 1 x R (these are 
sometimes called Steiner circles) where each ixR lies in the ideal boundary of a hyperbolic 
halfplane bounded by a, two such lines x x R, i/xR are at distance 9 in the S 1 factor 
if they meet at da at angle 9. If R C S 2 — da, then define visualangle^ (R) = 9 G [0, 27r], 
where 9 = inf{6>2 — 9\ mod 2ir \ R C [9i, #2] x R}- The possible choice of or 2/T is made 
in the obvious manner. 

Proposition 2.3: Let Si,Sj be disjoint oriented geodesies in H 3 . Then there exists 
a smooth simple closed curve Kij in S 2 ^ separating dSi from dSj such that for k G 
{*,:?}, visualangle 5fc (k^) = C (distance Sj)). 

Proof: Let P be the orthocurve from Si to Sj. Consider the half-plane with boundary Si 
determined by Si and P, and the half-plane with boundary Sj determined by Sj and P. 
Allow these half-planes to expand into solid angles at the same rate. (A solid angle is a 
closed set in B 3 = H 3 U bounded by two hyperboic halfplanes which meet along a 
common geodesic.) At first, the four half-planes that bound these solid angles intersect in 
H 3 , but at some angle 9 these half-planes intersect only at infinity (that is, in S^). By 
reasons of symmetry they intersect in two or four points (four points of intersection occur 
when Im(distance(5i, Sj)) is tt/2 or 3n/2). 

Let Si, Sj be the solid angles which exist at angle 9. Let = Sk H S^ for k G {i,j}- 
Let be a simple closed curve in Tj n Tj which separates dSi from 85 j. By construction, 
for k G {z, j} visualangle^^^ (/Cy) = 9. See Figure 2.2. 

To complete the proof of the Proposition, we now show that 9 = C(distance(5i, Sj)). To 
do this, we use hyperbolic trigonometry on a degenerate right-angled hexagon in H 3 . Fol- 
lowing [F], a degenerate right-angled hexagon is a 5-tuple of oriented geodesies Si, ■ ■ ■ , ^5 
in H 3 such that Si is orthogonal to Si+\ and Si and S5 limit at a common point at 
infinity. These oriented geodesies give rise to complex numbers cr , cr 2 , 03, 04. o"o = if 
the axes Si and £5 either both point into Sq or both point out of Sq. Otherwise ao = ni. 
For k G {2, 3, 4} cr^ = e if an e-translation of the oriented geodesic Sk takes the ori- 
ented geodesic S^-i to the oriented geodesic S^+i ■ By [F; pg. 83] we have the following 
Hyperbolic Law of Cosines. 

cosh(cr ) = cosh(cr 2 ) cosh(<7 4 ) + sinh(<j2) sinh(a 4 ) cosh(cr 3 ). 

We work in the upper-half-space model of hyperbolic 3-space, and normalize so that 
the ortholine from Sj to Si is B (thus Si intersects B above Sj), while the oriented axis Si 
is C. B will be S3, while the oriented geodesies Si and Sj will be S2 and £4, respectively. 

Of course, u + iv = distance^, Sj). If — tv/2 < v < then the intersection points at 
infinity occur in the second quadrant and the fourth quadrant (see Figure 2.3a). For con- 
venience, we work with the point in the second quadrant and send (unique) perpendiculars 
from it to the geodesies Si and Sj. The perpendicular to Si will be oriented towards Si and 
then denoted Si, while the perpendicular to Sj will be oriented away from Sj and then 
denoted S5. The intersection point at infinity (in the second quadrant) is So. 

This is the proper set-up for applying the (degenerate) Hyperbolic Law of Cosines (see 
Figure 2.3b). Note that a 3 = —(u + iv), and ao = iTv. By symmetry 02 = 04 = (a + i(3)/2 
where (a + i/3)/2 is distance(Si, S3). Plugging into the Law of Cosines, using a half-angle 
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Here dist(5j , 6i) is approximately .9+2/5 
\ and C(dist(Sj , 5i)) is approximately 2/3 / 

\ Figure 2.2 / 

formula (cosh(2z) = 2 cosh 2 (z) Al = 2 sinh 2 (z) + 1), solving for cosh(a + i(3), and taking 
the Arccosh, we get the desired result. Note that the visual angle in this set-up is —f3, 
thus necessitating taking the absolute value. 

When < v < n/2 our 2 intersection p7nnts~oTcux in the first and third quadrants, 
and we carry out the same procedure. This time S2 and S4 are traversed in the direction 
opposite to their orientations (the attendant changes in sign drop out though). In this 
case, the visual angle is (3. 

The cases —n < v < —tt/2 and n/2 < v < n reduce to the previous cases after 
adding or subtracting n. The formula in Definition 2.1 is then obtained after noting that 
cosh(z ± in) = — cosh(z). I 

Definition 2.4: Let 5 be a simple closed geodesic in the closed orientable hyperbolic 
3-manifold N. Let {5;};> be the lifts of 5 to H 3 . For each Tri^AQ-orbit of unordered 
pairs (Si 7 5j) choose a representative where i = 0. If Redistance(5 , 5i) < ln(3)/2, then 
let KQj be a smooth simple closed curve in S 2 separating 88q from dSj such that for 
k G {0,j}, visualangle 5fc (^oj) = C (distance (5o, If Redistance(5o 7 ^i) > ln(3) /2, then 
let Koj be the Dirichlet insulater, i.e. the boundary of the geodesic plane orthogonally 
bisecting the orthocurve between 5o, Sj. 

In either case, extend the collection 7Ti(iV)-equivariantly to a family {Kij} defined for 
all i, j. This is the Corona family for 5. 
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Si — S 2 




Proof: We check that {/%} satisfies the various conditions of Definitions 0.4-0.5 of [G]. 

i) By construction, Kij separates ddi from 85 j and {fty} is 7Ti(iV)-equivariant. Because 
for k G Sk — visualangle(Kij) < tt, satisfies the convexity condition. 

Modulo the natural action of tti(N) on k^, there are only finitely many insulators Ky 
which are not Dirichlet insulators. Therefore, for fixed i, there exist only finitely many 
such that diam(Kij) > e. This establishes local finiteness. | 

ii) No trilinking follows immediately from ii). | 

Definition 2.6: If 6 is a simple closed geodesic in the hyperbolic 3-manifold N, define 
maxcorona(5) = max{C(5o, Sj) \ j > 0} 

Remark 2.7: a) It seems possible that the Dirichlet insulator family associated to the 
geodesic 5 G N may be non-coalescable, while the Corona insulator family is coalescable 
and conversely. 

b) If tuberadius(5) > ln(3)/2, then both the Dirichlet insulator and Corona insulator 
families coincide and by Lemma 5.9 of [G] they are non-coalescable. 
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Proposition 2.8: Let 5 be a shortest geodesic in the closed orientable hyperbolic 3- 
manifold N. Then either the Corona insulator family is noncoalescable or there exists a 
2-generator subgroup 67 of wi (N) with generators /, w such that the parameter associated 
to {G, f, w} lies in the sub-box Xq = Xo a U Xq^ C W. 

Proof: Let 6 be a shortest geodesic in N, a closed orientable hyperbolic 3-manifold. By 
Corollary 1.29 and Lemma 2.5 either 6 has a non-coalescable insulator family or length(5) 
lies in C(Xk) for some k G {0, 1, . . . , 6} and maxcorona(5) > 27r/3. Let 67 be a 2-parameter 
subgroup of tti(N) generated by / and g, where / corresponds to 5 and has axis B and 
g maximizes C(g(B), B). We will show that if C(g(B),B) > 27r/3, then the parameter 
associated to {67, /, g} lies in the sub-box Xq. 

In fact we will show that if {6, /, g} is any torsion-free subgroup of Isom(H 3 ) where / 
is a length-minimizing loxodromic element with axis B and g maximizes C(g(B), B) where 
C(g(B),B) > 2tt/3, then the parameter /3 = (Lp 7 Dp,Rp) G exp _1 (X ). 

It follows as in the first paragraph that Lp e C(Xk) for some k G {0, 1, 2, 3, 4, 5, 6}. 
Also Dp is subject to nontrivial constraint. For example if Lp e C(X\), then Dp must lie 
in the decorated region of Figure 2.4, because Redistance((7(i?), B) > mind(Xfc) > 1.059 
where min^Xfc) is the minimal d value in Xf~. In fact the disjointness of the L(Xfc)'s 
implies that if (3 e S fl then tuberadius(5) > min d (X fc )/2, where 5 corresponds to the 
element /. Finally C(g(B),B) > 2tt/3 implies that Dp cannot lie in the decorated region 
of Figure 2.4. 



3 
2.5 
2.0 

a 

a 1.5 

1.0 
0.5 


Our proof is now similar to 8 the proof of Proposition 1.28*. We partition°the initial box 
W into sub-boxes Wi and eliminate Wj if any of theFfeMowing conditions hold (for clarity, 
the conditions have been translated into "pre-exponentiated" form). 

a) There exists no f3 G Wi such that length(/ / g) G C(X k ). 

b) Wi has some L values in C(Xk) but there exists no f3 G Wi such that distance 
{wp{B), B) G decorated region for that fc. 




21 



DAVID GABAI, G. ROBERT MEYERHOFF, AND NATHANIEL THURSTON 



c) There exists a killerword h in /, w, such that Relength^) < Relength^) 
and hp ^ id for all (3 GWj. 

d) There exists a killer word /i in /, w, f -1 , w~ l such that 

C(distance(/i /3 (S),S)) > C (distance (^(.B), S)) 

and hp{B) ^ S for all jSeWj. 

We have two files that contain the decomposition of W into sub-boxes and associated 
conditions/ (killer words). The program fudging checks that these files do indeed eliminate 
all of W — Xq. Fudging analyzes the cases k = 1, . . . , 6 all at once. | 

Remarks 2.9: i) Note that in practice the proof of Proposition 2.8 requires working in a 
considerably smaller parameter space than that of Proposition 1.28. Condition a) implies 
that the parameter space is "(2 + e)-complex dimensional" and condition b) implies that 
one of these parameters is greatly constrained. This suggests why it took so much longer 
to come up with the partition and the associated killer words for Proposition 1.28. In fact, 
it took roughly 1500 CPU days to find the partition and the associated killer words for 
Proposition 1.28, versus roughly 2 CPU days for Proposition 2.8. Here, the term "CPU 
day" refers to 24 hours of running an SGI Indigo 2 workstation with an R4400 chip, and 
the estimate of 1500 CPU days refers to 15 to 20 such machines running 80 to 90 percent 
of the time for 3 to 4 months. 

ii) We took pains to make fudging as similar to verify as we could, thereby lessening 
the amount of analysis needed to show the veracity of fudging. Appendix 2 contains those 
sections of fudging that differ from corresponding sections of verify. 

iii) When working with exponentiated co-ordinates (that is, in W rather than V) the 
Corona function changes as follows. Let X = exp(ct + ij3) and U = exp(u + iv), then the 
Corona function formula 

4 

cosh(a + i(3) = 1 — 

1 + cosh(w + iv) 

becomes 

X + X' 1 _ j 4 

2 ~ 1 ~ 1 + (U + U~ 1 )/2 

It is a pleasant exercise to solve this, and we find that 

(U 2 -6U + 1)±4(U-1)V=U 
(U + l) 2 

the two answers so gotten are reciprocals, which means their associated arguments are 
opposites. In fudging, the exponentiated version of the Corona function is the function 
horizon (ortho), which takes in U =ortho and computes the associated X value. (3, the 
argument of X, is implicitly gotten in the function larger-angle. 

iv) It is possible that by working purely in the context of the Corona function, rather 
than first working with Redistance and attempting to prove Proposition 1.28, the computer 
proof can be simplified. We started this project with the naive idea that perhaps Vol3 was 
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the only manifold whose shortest geodesic did not have a ln(3)/2 tube. The remarkable 
fact that this naive idea is almost correct accounts for the fact that a proof of Theorem 
0.2 can be obtained with only the mild extra effort detailed in this chapter and the next. 
For < v < tv/2 C(u, v ) = (3 where (3 satisfies 

A 

cosh(a + i(3) = 1 



1 + cosh(w + iv) 
For — 7r/2 < v < C(u, v) = —(3 where (3 satisfies 

4 

cosh(a + i(3) = 1 -. r. 

1 + cosh(w + iv) 

For 7r/2 < v < 7r C(u 7 v) = —f3 where (3 satisfies 

4 



cosh(a + i(3) = 1 



1 + cosh(tt + i(v — 7r)) ' 
For — iv < v < —iv/2 v) = (3 where (3 satisfies 

4 

coshfai + id) = 1 ; ; — . 

1 + cosh(-u + i(v + 7r)) 

In all of these cases, we normalize the sign of (3 by requiring that a be positive. 

Lemma 2.10: i) if h = wFwfwfWfwf then Relength(/i a ) < Relength(/ a ) for all pa- 
rameters a in exp _1 (Xo) 

ii) if /12 = wFwfwwfwFw then Relength(/i a ) < Relength(/ a ) for all parameters a 
in exp _1 (Xo) 

iii) \R\ < \L/2\ throughout exp" 1 ^). 

Proof: This easy calculation was carried out by verify. I 



Chapter 3: Vol3 

The two main results of this chapter are 

Proposition 3.1: If iV is a closed orientable hyperbolic 3- manifold, then either iV is Vol3, 
the third smallest known closed orientable hyperbolic 3-manifold, or maxcorona(5) < 27r/3 
for 6 a shortest geodesic in N. 

Proposition 3.2: Any shortest geodesic in Vol3 satisfies the insulator condition. 

Remark 3.3: Vol3 is the third smallest known closed orientable hyperbolic 3-manifold. 
Topologically Vol3 is (3,1) surgery on manifold m007 in the census of cusped hyperbolic 
3-manifolds (see [W]). It is also (-3,2) (-6,1) surgery on the "left-handed Whitehead link", 
link 5| in the standard knot tables. Snappea gives an experimental proof that Vol3 is 
hyperbolic and that its volume is that of the regular ideal 3-simplex. A rigorous proof can 
be found in [JR]. Previously, Hodgson- Weeks [HW1] had found an exact Dirichlet domain 
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for Vol3, that is, the face pairings were expressible as explicit matrices with coefficients 
in a finite extension F of Q and they obtained equations in F for the various faces. See 
Remark 3.14. 

Remarks 3.4: i) Idea of Proof of Proposition 3.1: If N has maxcorona(5) > 2/t/3, then 
it must have an (L,D,R) parameter in the region 1Z = exp _1 (Xo), (Xq is defined in 
Proposition 1.28). A geometric argument (Lemma 3.7) which utilizes Lemmas 3.5 and 3.6 
shows that R = and an algebraic argument (Lemmas 3.8, 3.9, 3.11, 3.12) shows that 
L = D = uj, where exp(u;) is a root of the polynomial 1 + 2d + 6d 2 + 2d 3 + <i 4 . This implies 
([HW1] or [JR]) that g = tti(Vo13) and hence N is covered by Vol3. By [JR], N = Vol3. 

ii) The proof of Proposition 3.2 will be somewhat similar to the proofs of Propositions 
1.28, 2.8. 

Lemma 3.5: Let K = exp _1 (X ) If a = (L, D, R) = (I + it, d + ib,r + ia) G K n T, then 
f a ,w a satisfy the relations 

i) wFwfwfWfwf 

ii) wFwfwwfwFw 

Proof: In i), ii) above and what follows below we supress the subscripts a. Also W 
(resp. F) denotes w -1 (resp. / _1 ). because i), ii) are cyclic permutations of the quasi- 
relators r2, t\ corresponding to the Xq sub-box of Proposition 1.28, it follows that if 
h = wFwfwfWfwf or wFwfwwfwFw, then Relength(/i) < Relength(/) throughout 7Z. 
Since a G T, / is a shortest element and so h = id. I 

Lemma 3.6: The following substitutions give rise to three sets of new relators. 

a) In i), ii) replace / by w and replace w by /. 

b) In i), ii) replace / by F. 

c) In i), ii) replace w by W. 

Proof: Note that replacing / by w necessitates replacing F by W, and so on. 

a) First, a cyclic permutation of relator i) gives relator i) with / replaced by w and 
w replaced by /. Second, one readily obtains the relator fW fwf fwfW f from i), ii), 
because fWfwf = (wFwfw) -1 = wfwFw = (fwfWf) -1 where the first and third 
equalities follow from i) and the second from ii). 

b) Again it is routine to obtain b) from i), ii). 

c) Conclusion c) follows from a) and b). I 

Lemma 3.7: If (L, D,R) e1Z(~)T, then R = 0. 

Proof: With respect to the parameter (L, D, R) let G be the group generated by /, w. 
Figure 3.1 shows a schematic picture of geodesies B, W(B), w{B), f(w(B)), f(W(B)). Also, 
it shows the images of the orthocurve O from W(B) to B after translation by w, f and fw. 
Finally, it shows the orthocurves 0\ from fW(B) to W(B) and O2 from w(B) to fw(B). 
Note that Figure 3.1 displays the situation where Re(i?) > 0. It is also apriori possible 
that O2 might intersect w(B) on the positive side of w(0). There are other similar possible 
inaccuracies. 

ui = fwffw G G sends geodesic W(B) to fw(B) and 02 = wFFwF G G sends the 
geodesic fW(B) to w(B). Now a^Vi = fWffWfwffw is a relator of G, since it is a 
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cyclic permutation of relator ii) with / replaced by w and w replaced by /. Thus o\ = 02 
and hence <j\(0\) = O2. 

Figure 3.1 

As in Fenchel [F] a right angled hexagon consists of a cyclically ordered 6-tuple of 
oriented geodesies in H 3 such that A^ intersects Ai+i (mod 6) orthogonally. 

Each "edge" of the hexagon is labelled by the complex number where one obtains A^+i 
from A;_i by translation of H 3 along the oriented geodesic A^. 

(3.1) The effect of reversing the orientation of A^ is to change to — e^, e^+i to ei+i+7rz, 
and ei-i to e^-i + ivi. 

Figure 3.1 gives rise to the two right-angled hexagons drawn in Figure 3.2. (Figure 3.2a 
may be inaccurate for the following reason. It is not clear whether the head of Oi should be 
placed in front of the tail of O or behind the tail of O. A similar statement holds for the tail 
of Oi and for Figure 3.2b.) Assume that in Figure 3.2a, Ai corresponds to B and the edges 
are cyclically ordered counterclockwise. Then = D, e\ = L, = —D. We now show 
that if es has value c, then has value c + ni. Observe that there is an order-2 rotation 
r of H 3 about an axis orthogonal to B which reverses the orientation on B and takes 
the oriented orthocurve O to the oriented orthocurve f(0). Since distance(VF(S), B) = 
distance(/W(S),S) it follows that r(fW(B)) = -W(B) and r(W(B)) = -fW(B) where 
the - sign indicates that the orientation has been reversed. This in turn imples that 
t(Oi) = — 0\ and therefore using (3.1) that = e$ + ni. 

Let (j) be the isometry of H 3 which is an r + i(n + a) translation of B. Thus 4>(B) = 
B, (f)(0) = -w(0), </>(f(0)) = -fw(0). This implies that <1>{W{B)) = w(B) and 
4>(fW(B)) = fw(B) which in turn implies that 4>(Oi) = — O2. If the hexagon of Fig- 
ure 3.2b with edges A' l7 • • •, Ag is counterclockwise cyclically oriented so that X[ denotes 



25 



DAVID GABAI, G. ROBERT MEYERHOFF, AND NATHANIEL THURSTON 




Figure 3.2 
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the oriented geodesic B, then again using (3.1) it follows that 



C+7I1 



O* 



r+ia 
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c + ni. 
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tr 

on Figure 3.3. Here O* is G translation of Oj with an endpoint on B. Two such translates 
appear in Figure 3.3, one where O* points iJfftguig 3aiid one where O* points out. Call the 
former (resp. latter) the pointing in (resp. out) O*. Actually Figure 3.3 includes 3 more 
relations. Because the oriented 0\ is a G— translate of the oriented 2 and / is a primitive 
element of G which fixes B 1 it follows that 

distance( (pointing in OjJ"), (pointing in O^)) = (mod L) 

and 

distance( (pointing out O*), (pointing out 0%)) =0 (mod L). 
Finally distance^ (O), O) = R. 
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We therefore obtain the following two equations 

(3.2) c - R + L + c + ni = (mod L) c-L + R + c + ni = (mod L) 
and hence 

(3.3) 212 = (mod L) Ac = (mod L). 

Using the V and R' ranges for Xq provided in Proposition 1.28, it is easy to compute that 
for each element of K = exp _1 (X ), |Re(i?)| < |Re(L)/2|. It then follows that R = for 
(L, D, R) G K n T. I 

Lemma 3.8: If (L, D,R) eTZdT, then L = D. 

Proof: We will now use the exponential coordinates / = exp(L), d = exp(D). These /, d 
should not be confused with the I + it,d + ib used above. In the following calculations 
Mathematica [Math] was used to perform matrix multiplication of 2 x 2 matrices with coe- 
ficients rational functions in the variables y/l, \fd. We will follow Mathematical notation 
(for example, in Mathematical notation the 2x2 matrix (a^) is {{an, a^}, {021, ^22}})- 
In particular, plugging in R = in Lemma 1.24 we get the following matrix representations 
of /, F, w, W. Because the R term drops out, we can express the matrices of w and W as 
functions of d alone. 

/[/_] = {{Sqrt[Z],0},{0,l/Sqrt[Z]}} 
F[L] = {{l/Sqrt[/],0},{0,Sqrt[/]}} 

w[d_] = {{l/2(Sqrt[d] + 1/Sqrt[d]), l/2(Sqrt[d] - 1/Sqrt[d])}, 
{l/2(Sqrt[d] - 1/Sqrt[d]), l/2(Sqrt[d] + 1/Sqrt[d])}} 

W[dJ\ = {{l/2(Sqrt[d] + 1/Sqrt[d]), l/2(-Sqrt[d] + 1/Sqrt[d])}, 
{l/2(-Sqrt[d] + 1/Sqrt[d]), l/2(Sqrt[d] + 1/Sqrt[d])}} 

Let Y be the relator F.w.f.w.f.W.f.w.f.w, which is a cyclic permutation of relator i) 
of Lemma 3.4. Multiplying this product of 10 matrices we obtain the following matrix for 
Y which we know is I. 

Y[L, d_] = {{((1 + d) * (1 - 2 * d 2 + d 4 + 8 * d * I - 16 * d 2 * I 
+8 * d 3 * Z - 2 * Z 2 + 4 * d * Z 2 - 4 * d 2 * Z 2 + 4 * d 3 * Z 2 - 
2 * d 4 * I 2 + Z 4 + 4 * d * Z 4 + 6 * d 2 * Z 4 + 4 * d 3 * Z 4 
+d 4 *Z 4 ))/(32*d( 5 / 2 )*Z( 3 / 2 )), 

((-1 + d) * (1 + 4 * d + 6 * d 2 + 4 * d 3 + d 4 + 4 * d * Z+ 
8 * d 2 * Z + 4 * d 3 * Z - 2 * Z 2 - 12 * d 2 * Z 2 - 2 * d 4 * Z 2 + 

4 * d * Z 3 + 8 * d 2 * Z 3 + 4 * d 3 * Z 3 + Z 4 + 4 * d * Z 4 + 
6 * d 2 * Z 4 + 4 * d 3 * Z 4 + d 4 * Z 4 ))/(32 * d { W * Z^ 2 ))}, 
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{((-1 + d) * (1 + 4 * d + 6 * d 2 + 4 * d 3 + d 4 + 4 * d * 1+ 
8 * d 2 * I + 4 * d 3 * I - 2 * I 2 - 12 * d 2 * I 2 - 2 * d 4 * l 2 + 

4 * d * I 3 + 8 * d 2 * I 3 + 4 * d 3 * I 3 + I 4 + 4 * d * Z 4 + 
6 * d 2 * I 4 + 4 * d 3 * I 4 + d 4 * / 4 ))/(32 * d^'V * /( 3 / 2) ), 



((1 + d) * (1 + 4 * d + 6 * d 2 + 4 * d 3 + d 4 - 2 * I 2 + 4 * d * l 2 - 
4 * d 2 * I 2 + 4 * d 3 * I 2 - 2 * d 4 * I 2 + 8 * d * l 3 - 
16 * d 2 * I 3 + 8 * d 3 * I 3 + I 4 - 2 * d 2 * I 4 + d 4 * I 4 )) 
/(32*rf( 5 / 2 U/( 3 / 2 ))}} 

Since G is generated by /, an L translation along B, and w, a D translation along 
C, it follows from Lemma 3.6 that the relation Y = I holds with / and d switched. Thus 
= Y12, and = Y 12 (with /, d switched) which implies that 

= (l + 4*d + 6*d 2 +4*d 3 + d 4 + 4*d*Z + 8*d 2 *Z+ 
4 * d 3 * I - 2 * I 2 - 12 * d 2 * I 2 - 2 * d 4 * l 2 + 
4*d*l 3 + 8*d 2 *l 3 + 4*d 3 *l 3 + l 4 + 4*d*l 4 + 

6 * d 2 * I 4 + 4* d 3 * I 4 + d 4 * l 4 )- 
(l + 4*l + 6*l 2 + 4*l 3 + l 4 + 4*l*d + 8*l 2 *d+ 

4 * I 3 * d - 2 * d 2 - 12 * I 2 * d 2 - 2 * I 4 * d 2 + 
4 * I * d 3 + 8 * I 2 * d 3 + 4 * I 3 * d 3 + d 4 + 4 * I * d 4 + 
6 * I 2 * d 4 + 4 * I 3 * d 4 + I 4 * d 4 ) = 

4 * (1 + d) 2 * (1 + 1) 2 * (-d + 1) * (-1 + d * 0- 

This implies that d = I and hence D = L or we obtain one of the following solutions 
which contradicts the condition Re(D) > 0. The solution d = — 1 implies D = ln(d) = 
ln(— 1) = 7ii. The solution d = 1 implies D = 0. The solution d*l = 1 implies that d = l/l 
and hence D = ln(d) = ln(l/Z) = - ln(Z) = — L and hence Re(L>) < 0. | 

Lemma 3.9: If (L, D,R) elZr\T, then <i = exp(D) is a root of the polynomial 

l + 2*d + 6*d 2 + 2*d 3 + d 4 . 

Proof: The equation Yi 2 = yields 
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= l + 4*d + 6*d 2 +4*d 3 + d 4 + 4*d*l + 8*d 2 *l + 4*d 3 *l-2*l 2 - 
12 * d 2 * I 2 - 2 * d 4 * I 2 + 4 * d * I s + 8 * d 2 * Z 3 + 4 * d 3 * Z 3 + I 4 + 4 * d * Z 4 + 

6 * d 2 * Z 4 + 4 * d 3 * Z 4 + d 4 * Z 4 . 

Setting Z = d we obtain 

= 

l+4*d+8* d 2 + 12 * d 3 -2* d 4 + 12 * d 5 + 8* d 6 + 4* d 7 + d 8 = 
(I + 2 * d - 2 * d 2 + 2 * d 3 + d 4 ) * (I + 2 * d + 6 * d 2 + 2 * d 3 + d 4 ). 

On the other hand setting Z = d in the equation Yn = 1 we obtain 

32d 5 = (1 + d)(l + 4 * d 2 - 12 * d 3 + 6 * d 4 + 8 * d 5 + 4 * d 6 + 4 * d 7 + d 8 )] 

and so 

= (-1 + d) * (1 + 2 * d + 6 * d 2 + 2 * d 3 + d 4 ) * (-1 + 4 * d 3 + d 4 ) 
The only common solutions to these equations are the roots of the equation 
(3.4) (l + 2*d + 6*d 2 + 2*d 3 + d 4 ) | 

Remark 3.10: The two facts Y = I, and Z, d can be switched in Y follow from relator i) 
of Lemma 3.4 and R = 0. Relator ii) was used in the proof that R = 0. 

Lemma 3.11: The roots of 1 + 2 * d + 6 * d 2 + 2 * d 3 + d 4 are 

(-1 + i* V3)/2 -(-6-2*i* v / 3) (1/2) /2 

(-1 + i* VS)/2 + (-6 - 2 * i * V3) {1/2) /2 
V3)/2 - (-6 + 2 * i * y/3) {1/2) /2 
(-1 - i * >/3)/2 + (-6 + 2 * i * v / 3) (1/2) /2 

I 

Remark 3.12: Note that if x is a root of 1 + 2 * cZ + 6 * d 2 + 2 * d 3 + cZ 4 , then so are x, 1/x 
and 1/x. 

Lemma 3.13: If (L, D,R) gTZDT, then 

D = L = ln((-l + i * V3)/2 - (-6 - 2 * i * >/3) (1/2) /2) 
= 0.83144294552931 - 1.945530759503636 * i. 

Proof: The other 3 solutions are — L, Re(L) — Im(L), —Re(L) + Im(L) and lie outside 1Z. 
The above solution lies in 1Z. I 
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Remark 3.14: In our language, Hodgson and Weeks [HW1] knew that 7Ti(Vol3) was 
generated by /, w with D = L = uj and R = 0. Also that the various solutions of (3.4) 
corresponded to symmetries of Vol3. 

Corollary 3.15: Vol3 is the unique hyperbolic 3-manifold with ^associated parameter 
values in S fl X . 

Proof of Proposition 3.1: By Proposition 2.8 if maxcorona(5) > ln(3)/2, then the 
parameter for a 2-generator subgroup G of tti(N) lies in S HXq. Since G = 7Ti(Vol3), N is 
covered by V613. By Reid (see [JR]) Vol3 only covers Vol3. Therefore N =Vol3. I 

Proof Proposition 3.2: We show that the Dirichlet insulator associated to a shortest 
geodesic S in N satisfies the insulator condition if C(0(2)) < 113.16 degrees. Here 0(2) = 
distance(£?, second closest lift of S) and 0(1) = oj corresponds to the closest lift. This 
actually requires that we be a bit more precise, so we insert the following definition. 

Definition 3.16: Let S be a geodesic in the hyperbolic 3-manifold N and {Si}i> be the 
preimages of S in H 3 . Define an equivalence relation on {di}i>i by saying that Si = Sj if 
either 

i) there exists a g G tti(N) such that g(So) = Sq and g(Si) = Sj 

or 

ii) there exists g G tti(N) such that g(5i) = S and g(So) = Sj. 

Call an equivalence class an orthoequivalence class. Consider the collection O of com- 
plex numbers {distance^, $i) | * > and only one Si is represented in each equivalence 
class } Now order O to obtain the ortholength spectrum of S, {0(1), 0(2), • • •} where % < j 
implies Re(0(z)) < Re(0(j))- 

The based ortholength spectrum consists of all distinct pairs of complex numbers of 
the form {(distance(C, ortholine from B to g(B)),dista,nce(B,g(B)) | g G n 1 (N)}. 

The geodesic plane midway between B and w~ 1 (B) intersects the B — C plane in a 
geodesic E at distance Re(u;)/2 = Re(length/)/2 from B. Conjugate tci(N) so that / is 
still an uj translation of B and so that one endpoint of E is at (1,0) and the other endpoint 
is at (x,0), where x > 1. By [Beardon; p. 166 ] x < 23.815. Thus the Dirichlet insulator 
A (resp. w(\) between B, w~ 1 (B) (resp. B, w(B)) is symmetric about the and 
lies within the circle passing through (1,0), (23.815,0) (resp. (-1,0), (-23.815, 0)). By [G; 
Lemma 4.7] this circle takes up a visual angle of less than 133.68 degrees. Now / is the 
composition of an exp(Re(o;)) homothety centered about the origin and an Im(u) radian 
w -111.4707 degree rotation. Because exp(4 * Re (a;)) > 27.82 > 23.815 it follows that 
the geodesic E is taken "beyond" E by w 4 . In fact, we see that A fl (f n (X) U f n w(X)) = 
if \n\ > 4. Therefore if there was a trilinking among three insulators associated to the 
orthoclass of w(B), there would be a tri-linking involoving 3 circles from the collection 
{f n w(X),f n (X) | -3 < n < 0}. 

Since is rotated by at most 111.48 degrees under / ±1 , and A takes up less than 
133.68 degrees B- visual angle it follows that / _1 (A) U A take up less than 133.68 + 111.48 = 
245.06 degrees B— visual angle. Similar arguments show that if A and one of f n (X) or 
f m w(X) nontrivially intersect, then the union cannot take up more S-visual angle, in fact 
except for / ± (A), it takes up less. See Figure 3.4 which shows the union of {f n w(X), f n (X) | 



30 



HOMOTOPY HYPERBOLIC 3-MANIFOLDS ARE HYPERBOLIC 



— 3 < n < 0}. Therefore, if the union of three such circles was connected, they would take 
up at most 133.68 + 2(111.48) = 356.64 degrees of S-visual angle and hence would not 
create a trilinking. 




Now suppose that C(0(2)) < 113.16. Thus an^. insu|ator associated to a translate of B 
not in the 0(1) class would take less than 113.16 degrees visual angle. There are three cases 
to consider: trilinking involving no 0(1) insulators; trilinking involving exactly one O(l) 
insulator; trilinking involving exactly two O(l) insulators. The case of no O(l) insulators 
cannot occur because 3(113.16) < 360. The case of exactly one O(l) insulator cannot 
occur because B-visual angle(0(l)) + 2(0(0(2)) < 133.68 + 2(113.16) = 360. Finally, the 
case of exactly two O(l) insulators cannot occur, because the two O(l) insulators would 
take up at most 245.06 degrees of visual angle (see the above paragraph) and adding to 
this the fact that 0(0(2)) < 113.16 produces less than 360 degrees. 

The proof of Proposition 3.2 is now completed by the following 

Lemma 3.17: If 5 is a shortest geodesic in Vol 3, then 0(0(2)) < 113.16. 

Proof: Step 1: Re(0(2)) > Re(0(l)) 

Proof of Step 1. If Re(0(2)) =Re(0(l)), then it follows by Proposition 3.1 that 
0(2) = O(l) = lo. Furthermore the collection of orthocurves must appear symmetrically 
on the geodesic B, i.e. if distance(0(2)— ortholine, C) = x+yi, with < x, x minimal, then 
there is an ortholine at distance n(x + yi) with n E Z along B from C whose corresponding 
ortholength is uj. Use the fact that if v E 7Ti(Vol3) is the element with distance(i> (B) : B) = 
0(2) as above, then the group generated by v, f is conjugate to the group generated by 
wj. 

If x = 0, and y > is minimal, then y = n/m. Using Fenchel's law of cosines [F; p. 
83 ], it follows that if m = 2, then there exists an uj ortholength thereby contradicting 
Proposition 3.1. See Figure 3.5. If m > 2, then a similar argument shows that there exists 
a real ortholength less than Re(w). 

If x > choose m minimal so that mx =Re(u;). By replacing y by y + ni, if necessary, 
we can assume that m(x + yi) = uj. Therefore an x + iy translation r along B descends 
to a Z/m action 4> on Vol3. Any lift of (p n is a conjugate of r n which is fixed-point free or 
the identity. This contradicts the fact that Vol3 only covers Vol3. I 
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Cosh(P)=Cosh 2 (co) + Sinh 2 (03)Cosh(7ii) 

Step 2: C(0(2)) < 113.16 Figure3 . 5 

Proof of Step 2. We obtain this result with computer assistance in a manner similar 
to the proofs of Propositions 1.28 and 2.8. Our parameter space W is the usual initial box. 
As before, a parameter gives rise to a group G with generators /, v. (Here we reserve the 
letter w to denote one of the generators of 7ri(Vol3).) We consider the parameters U C W 
such that /, v generate a group G where / is of minimal length, and maxcorona(5/) = 
C(v) > 113.16. If Step 2 was false then U ^ for it contains the 2-generator subgroup of 
7Ti(Vol3) generated by /, v where distance^ (£?), B) = 0(2). We now show that U = as 
follows. 

We partition W into regions Wi such that each Wi can be eliminated for one of the 
following reasons, 
z) W, = X . 

a) There exists no (3 G Wi such that length(/ ( g) G C(Xq). In particular, length(/a) 7^ u 
throughout Wi. 

b) Wi has some L values in C(X ) but there exists no f3 G Wi such that the real 
part of distance^ (5), B) is greater than the minimum d value for Xq. In particular, 
Re(distance(w j a(S), B)) < Re (a;) throughout Wj. 

c) There exists a killerword h in /, w, to -1 , such that ^ 7^ id and Relength(/i^) < 
Relength(/ j3 ) for all (3eW t . 

It turns out that a Corona condition is not needed — c) is sufficient to generate all the killer 
words we need! 

There are files containing the partition of W and the associated conditions/killerwords, 
and fudging verifies that they indeed work. Fudging actually takes care of the sub-boxes 
Xq, . . . , X 6 all at once, and the associated files reflect that fact. 

As noted earlier, fudging works with the exponentiated versions of the above condi- 
tions. I 

Here is an experimental "proof of Lemma 3.17. In [HW2] an algorithm is given to 
compute, with multiplicities, the length spectrum of a hyperbolic 3-manifold M, given a 
Dirichlet domain for M. Weeks has observed [Weeks, personal communication] that a very 
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similar argument gives an algorithm to compute the based ortholength spectrum. In fact 
an analogue to Proposition 1.6.2 [HW2], (with an analogous proof) is the following 

Lemma 3.18 (Weeks): Let M be a closed orientable 3-manifold having a Dirichlet 
domain V with basepoint x and with spine radius r. Let 5 be a geodesic of length / + it. 
To compute all the based ortholengths of real length less than or equal to A with basing 
less than or equal to 1/2 from some point on So ( a preimage of S) it suffices to find all 
translates gV satisfying p(x,gx) < 2r + 2Arccosh(cosh(//2) cosh(A/2)). 

Proof: As in [HW2] we can assume that p(So,x) < r. The 0-basing on 5 wm be given 
by the oriented perpendicular P from Sq to x. Figure 3.6 shows that if there is a translate 
5i = g(S ) based at distance < 1/2, at real distance < A from 5 , then p(g(x), x) < A + /+2r. 
As in [HW2] an application of the hyperbolic cosine law yields the better estimate of the 
Lemma. | 

y 

\ Am W 
/ g(x) 

A 



Provided one is given a Dirichlet domlan, [HW2] gives an efficient algorithm to find 
these g : s. Finally to each such g one computes the basing and distance from g(So) to 8q. 

The collection of ortholengths with basings < 1/2 contains at least 2 representatives for 
each ortholength class, thus the Weeks algorithm can be used to give lower bounds on the 
various 0(i)'s. Snappea computes a Dirichlet domain for Vol3 with spine radius < 0.68. 
Because / < 0.83145, taking A = 1.24 we obtain 2r + 2arccosh(cosh(//2) cosh(A/2)) < 
2.89. This algorithm has been implemented on an undistributed version of Snappea, and 
provided the following estimates. Note that Re(0(2)) > 1.24 is sufficient to guarantee that 
C{0{2)) < 113.16 degrees. 

0(1) w .83144- 1.94553/, 

0(2) w 1.3170 -ttI, 

0(3) = 0(4) ss 1.4197 + 1.0963/, 

0(5) w 1.9769 - 1.2995. 

These estimates were found using a "tiling" of radius 3.00 > 2.89, which is sufficient for 
"proving" that for Vol3, Re(0(2)) > 1.24. | 

Remark 3.19: This experimental proof should be easily rigorizable by implementing 
Weeks's algorithm using exact arithmetic. 
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Chapter 4: Applications 

We provide some applications. First, we give a partial answer to Question 1.16. Second 
a lower bound for the volume of hyperbolic 3-manifolds is produced. Finally, we give a 
relationship between isotopic closed curves in a hyperbolic 3-manifold and essential links 
in B 3 . 

Theorem 4.1: If 5 is a shortest geodesic in a closed orientable hyperbolic 3-manifold N, 
then either tuberadius(5) > 1.059191579962 . . ./2 or N =Vol3. 

Proof: If tuberadius(5) = d/2 < ln(3)/2, then by Corollary 1.29 ii, d E 

[0.831426508686 . . . , 0.831461989726 . . .] U [1.068029907104 . . . , 1.068134862048 . . .]U 
[1.072078514724 . . . , 1.072622271888 . . .] U [1.094788113123 . . . , 1.095007375882 . . .]U 
[1.094973832380 . . . , 1.095231885426 . . .] U [1.059191579962 . . . , 1.060372139694 . . .]. 

If d < 1.059191579962..., then by Proposition 1.28 the parameter (L,D,R) e V asso- 
ciated to N must be in 1Z fl T = exp _1 (X H S). It then follows by Corollary 3.15 that 
N = Vol3. I 

Conjecture 4.2: If 6 is a shortest geodesic in a closed orientable hyperbolic 3-manifold 
N, then either tuberadius(5) > ln(3)/2 or N is one of six exceptional manifolds. 

Remark 4.3: Vol3 is one of these manifolds. As in Remark 1.32, three others are conjec- 
turally, s479(-3, 1), s778(-3, 1), u2678(2, 1). These correspond to sub-boxes X 5 ,X 2 ,Xi 
respectively. 

The conjectured fundamental groups of the six manifolds are < /, w; ri(Xk),r2(Xf~) > 
for k e 0, . . . , 6. (The groups for k = 5, 6 are isomorphic). Using [Be] one can get explicit 
Heegaard genus 2 descriptions of all the conjectured manifolds. 

One could prove Conjecture 4.2 by first showing that for each k, T fl Xk is a point 
Tfc. Let Mk be the hyperbolic 3-manifold which corresponds to Tfc. Second show that Mk 
nontrivially covers no 3-manifold. This procedure was carried out in Chapter 3 for the 
subbox Xq. 

Remarks 4.4: The previous best lower bound for the volume of hyperbolic 3-manifolds 
was on the order of 0.001 (see [GM1] and [M2]). Using the results of this paper and the 
method of [Ml] it is easy to improve this to 0.1. However, Gehring and Martin provide an 
improved tube-volume formula in [GM2] and we use their formula to get a lower bound 

of 0.16668 The Gehring-Martin tube-volume formula for manifolds (as opposed to 

orbifolds) is 

V(t) = v / 3tanh(t)cosh(2t)Arcsinh 2 (sinh(t)/cosh(2t)) 

where t is the radius of the embedded solid tube. Note that the length of the core geodesic 
is irrelevant. 

Theorem 4.5: ^-^Arcsinh^^p) = 0.16668 ... is a lower bound for the volume of closed 
hyperbolic 3-manifolds. 
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Proof: [GM2; Corollary 1.7] applies the tube- volume formula V(t) to tubes of radius at 
least ln(3)/2 and produces 

Vol(iV) > V(m(3)/2) = -^Arcsinh 2 (^?) = 0.16668. . . 

2v 3 5 

Corollaries 1.29 and 3.15 imply that if the tube radius of a shortest geodesic in N is 
less than ln(3)/2, then either N = Vol3 or N contains an embedded open tube of radius 
d/2, where d G Re(ln(D(Xfc))), about a geodesic of length / G Re(ln(£(X/ c ))) for some 
k G 1, . . . , 6. As in the Theorem 4.1 we compute that / > 1.059536368901 . . . 

In the first case, Vol(iV) = Vol(Vol3) = 1.01 ... , while in the second case I and d are 
bounded as follows: / > 1.059536368901... and d > 1.059191579962.... Plugging these 
into the tube-volume formula nl sinh 2 (<i/2) we get Vol(A^) > 1.02419 I 

Theorem 4.6: Let fei, ki be simple closed curves in N such that k\ is a geodesic. Then k\ 
is isotopic to ki if and only if as £? 3 -links q~ 1 (ki) is equivalent to q~ l (k,2) where q : H 3 — > N 
is the universal covering projection. 

Proof: Apply Corollary 5.6 of [G]. Recall that an equivalence between S 3 -links T, A is a 
homeomorphism of B 3 which takes V to A and fixes S 2 pointwise. I 

Remark 4.7: A similar argument extends Theorem 4.6 to homotopy essential links which 
lift to trivial £? 3 -links. The general case of Conjecture 5.5A of [G] is still open. 

Chapter 5: Conditions and Sub-Boxes 

As described in Section 1, our goal is to kill off all points in S C C 3 , but for com- 
putational reasons we will actually work with the box B D S. We will decompose B into 
a collection of sub-boxes such that each sub-box has an associated "condition" that will 
describe how to kill off that entire sub-box. The set-up for describing these sub-boxes will 
be given in Section 3. 

We now list the conditions used to kill off the sub-boxes. There are two types of 
conditions: the trivial and the interesting. The trivial conditions kill off sub-boxes in W 
by noting that the sub-box in question misses S. The interesting conditions are where the 
real work is done, and they require a killer word in /, w, w -1 to work their magic. 

To be consistent with the computer program verify we use the following notation. 
L' = z + iza, D' = zi + iz4, and R' = z 2 + iz 5 . Here (Z/, D' , R') G W and V = exp(L) = 
exp(/ + it), D' = exp(D) = exp(<i + ib), R' = exp(R) = exp(r + id). 

The Trivial Conditions 5.1: 

Condition 's' (short): Tests that all points in the sub-box have \z + iz$ | < 1.10274. 
This ensures that exp(Z) = |exp(L)| = \L'\ = \z + iz 3 \ < 1.10274 < exp(0.0978), and 
Definition 1.14 tells us that we are outside of S = exp(T) 

Condition '1' (long): Tests that all points in the sub-box have \z + iz 3 \ > 3.63201. 
This ensures that exp(Z) = |exp(L)| = \L'\ = \z + iz 3 \ > 3.63201 > exp(l. 2897845) and 
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we are outside of S. Here we have improved S to allow for the value computed in the proof 
of Proposition 1.13, rather than the cruder value 1.29. 

Condition 'n' (near): Tests that all points in the sub-box have \z± + iz 4 \ < 1. This 
ensures that exp(cf) = | exp(D)| = \D'\ = \z\ + iz 4 \ < 1 = exp(O) and we are outside of S. 
Actually, we used a stronger condition in Definition 1.14 (1/2 < d). 

Condition 'f (far): Tests that all points in the sub-box have \z\ + iz 4 \ > 3. This 
ensures that exp(<i) = | exp(D)| = \D'\ = \zi + iz 4 \ > 3 = exp(ln3) and we are outside of 
S. 

Condition 'w' (whirle big): Tests that all points in the sub-box have \z2 + iz§\ 2 > 
\zq + izs\ This ensures that exp(r) = |exp(i?)| = \R'\ = \z<i + iz§\ > a/I^o + iz%\ = 
a/ exp(/) = exp(7/2) and we are outside of S. In the computer program, this test requires 
a round-off error analysis. As such, the formula in the program is more complicated than 
expected. See Chapters 7 and 8 for a discussion of round-off error. 

Condition 'W (whirle small): Tests that all points in the sub-box have | Z2+iz^\ < 
1. This ensures that exp(r) = | exp(i?)| = \R'\ = \z2 + izs\ < 1 = exp(O) and we are outside 
of S. 

The Interesting Conditions 5.2: 

Condition 'L': This condition comes equipped with a killer word k in / and w, 
and tests that all points in the sub-box have | exp (length (k)) | < \L'\ = |exp(L)|, where 
length(/c) means the length of the isometry determined by k. This, of course, contradicts 
the fact that L is the length of the shortest geodesic. 

It is easy to carry out the test | exp(length(/c))| < \L'\ because Lemma 1.25a) can be 
used. Note that in verify the function which computes exp(length) is called length. 

Of course, Condition 'L' also checks that the isometry corresponding to the word k is 
not the identity. 

Condition 'O': This condition comes equipped with a killerword k in / and w, and 
tests that all points in the sub-box have | exp(distance(/c(S), B))\ < \D'\ = \ exp(D)|. This, 
of course, contradicts the "nearest" condition. 

It is easy to carry out the test | exp(distance(/c(£?), B))\ < \D'\ because Lemma 1.25b) 
can be used. Note that in verify the function which computes exp(distance(/c(S), S))is 
called orthodist. 

Also, Condition 'O' checks that the isometry corresponding to the word k does not 
take the axis of / to itself, bigskip Condition '2': This is just the 'L' condition without 
the "not-the-identity" check, but with the additional proviso that k is of the form f p w q . 
This ensures that k is not the identity, because for k to be the identity / and w would 
have to have the same axis, which contradicts the fact that d can be taken to be greater 
than or equal to 1/4. 

Condition 'conjugate': There is one other condition that is used to eliminate points 
in W. Following Definition 1.14 (and Lemma 1.15) we eliminate all boxes with < t < tt. 
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Of course, after exponentiating L = I + it, this corresponds to eliminating all boxes with 
zs > 0. Specifically, we toss all sub-boxes of W whose fourth entry is a 1. This condition 
does not appear in verify and fudging because it is applied "outside" of these programs, as 
described at the end of this Chapter. 

Construction 5.3: We now give the method for describing the roughly 930 million sub- 
boxes that the initial box W is sub-divided into. 

All sub-boxes are gotten by sub-division of a previous sub-box along a real hyper-plane 
mid-way between parallel faces of the sub-box before sub-division. Of course, these midway 
planes are of the form Xi = a constant. We use 0's and l's to describe which half of a sub- 
divided sub-box to take (0 corresponds to lesser Xi values). For example, describes the 
box Wn{(xo, xi, X2, £3, £4, £5) : x$ < 0}, 010 describes the box Wn{(£o; xi, X2, £3, £4, £5) : 
xo < 0, x\ > 0, X2 < 0}, and so on. 

In this way, we get a 1-to-l correspondence between strings and sub-boxes. If s is a 
string of 0's and l's, then let Z(s) denote the box corresponding to s. The range of values 
for the i-th coordinate in the sub-box Z(s) is related to the binary fraction 0.SiSi + Q . . . Sj + 6fc- 
The two sub-boxes gotten from subdividing Z(s) are Z(s0) and Z(sl). 

The directions of sub-division cycle among the various coordinate axes: the n-th sub- 
division is across the (n mod 6)-th axis. The dimensions of the top-level box W were 
chosen so that sub-division is always done across the longest dimension of the box, and so 
that all of the sub-boxes are similar. This explains the factor of 2^ 5- ^/ 6 in Definition 1.22. 

To kill a sub-box Z(s), the checker program has two (recursive) options: use a con- 
dition (and, if necessary, an associated killer word) to kill Z(s) directly, or first kill Z(s0) 
and then kill Z(sl). 

There is also a third option: don't kill Z(s), and instead mark s as omitted. Any 
omitted sub-boxes are checked with another instance of the checker program, unless the 
sub-box is one of the 7 exceptional sub-boxes (11 before joining abutters). 

Thus, a typical output from verify would be 

verifiedOOOOOOllllOlllllll - {0000001111011111110000000111101111111110} 

Which means that the sub-box Z(OOOOOOllllOlllllll) was killed except for its sub-boxes 
Z(0000001111011111110) and Z(OOOOOOllllOlllllllllO). The ouput 

verifiedOOOOOOllllOlllllllO- { }. 

and 

verifiedOOOOOOllllOlllllllllO- { }. 

shows that these boxes were subsequently killed as well, and thus the entire sub-box 
Z(OOOOOOllllOlllllll) has been killed. 

Instead of immediately working on killing the top-level box, we subdivide in the six 
co-ordinate directions to get the 64 sub-boxes 

Z(000000), Z(000001), Z(000010), Z(000011), . . . , Z(llllll). 
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We then throw out the ones with fourth co-ordinate equal to 1 (thanks to Definition 1.14 
and Lemma 1.15), leaving us with the 32 sub-boxes 

Z(000000), Z(000001), Z(000010), Z(000011), . . . , Z(111011). 

We then use verify to kill these, with the exception of the 7 exceptional boxes (11 before 
joining abutters together) listed in Proposition 1.28. 

The choices in verify are made for it by a sequence of integers given as input. The 
sequence of integers containing the directions for killing Z(000000) is contained in the file 
data6/000000. In such a sequence, tells verify to sub-divide the present box (by xi = c), 
to position itself on the "left-hand" box (xj < c) created by that sub-division, and to read 
in the next integer in the sequence. A positive integer n tells verify to kill the sub-box 
it is positioned at directly, using the condition (and killer word, if necessary) on line n in 
the 'conditionlist' file, and to then position itself at the "next" natural sub-box. —1 tells 
verify to omit the sub-box, and mark it as skipped (the sequence of integers used in killing 
the skipped box Z(s) is contained in a file data6/s).. 

The checker program verify, its inputs, and the list of conditions will be available from 
the Geometry Center. Details about how to get them can be found at 

http : //www . geom . umn . edu : /locate/HomotopyHyperbolic 

Similarly for the program fudging which is used on the 7 exceptional boxes. 

Example 5.4: To illustrate the checking in action, this is a (non-representative) example, 
which shows how the sub-box Z(s) (minus a hole) is killed, where 

s = 001000110001110111001111000101111111101111100111001111000001111011110111. 

The input associated with this sub-box is 

(0, 0, 0, 1929, 12304, 0, 0, 7, 0, 1965, 0, 1929, 1929, 1996, -1), 

which causes the program to kill Z(s) in the following fashion: 
kill Z(s): 
kill Z(s0): 
kill Z(s00): 

kill Z(s000) with condition 1929 = "L(FwFWFWfWFWFwFwfww)" 
kill Z(sOOl) with condition 12304 = "L ( FwfWFFWF wFwfwfW fwfw) " 
kill Z(s01): 
kill Z(s010): 

kill Z(s0100) with condition 7 = "L(w)" 

kill Z(sOlOl): 

kill Z(s01010) with condition 1965 = "L(fwFwFWFFWFwFwfwww)" 
kill Z(s01011): 
kill Z(sOlOllO) with condition 1929 
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kill Z(sOlOllO) with condition 1929 
kill Z(s011) with condition 1996 = "L(FwFwFWFWfWFWFwFwfww)" 
omit Z(sl) 
as shown in figure 5.1. 

Z(sl) is ignored, so the checker would indicate this omission in its report. In fact, 
Z(sl) is omitted entirely, since it is one of the 11 exceptional boxes. 

The use of condition "L(w)" so deep in the tree is unusual. In this case, it's because 
the manifold in the exceptional sub-box has length(f) = length(w), so that the program 
will frequently come to places where it can bound length(f) > length(w) nearby. 

One might wonder why the checker subdivides Z(s01011), since it's going to use the 
same condition to kill both halves. The reason is the error bound for Z(s01011) wasn't 
good enough to prove that the sub-box is killed directly. 

The binary numbers used by the computer require too much space to print. In the 
example calculation which follows, we instead use a decimal representation. Also, only 10 
decimals are printed, less accurate than the 53 binary digits used for the actual calculations. 

The sub-box Z(s01011) is the region where 

/ -1.381589027741 . . . < Re(L') < -1.379848991182 . . . \ 
-1.378124546093 . . . < Re(D') < -1.376574349753 . . . 

0.999893182771 . . . < Re(R') < 1.001274250703 . . . 
-2.535837191243 . . . < Im(L') < -2.534606799593 . . . 
2.535404997792 . . . < Im(D') < -2.534308843448 . . . 
V -0.001953125000 . . . < Im(R') < 0.000000000000 . . . / 



At this point, we would like to compute 

/, w, g = f~ 1 wf~ 1 w~ 1 f~ 1 w~ 1 fw~ 1 f~ 1 w~ 1 f~ 1 wf~ 1 wfww, length^), 

and so on. However, these items take on values over an entire sub-box and thus are 
computed via AffApprox's (first-order Taylor Approximations with remainder bounds), 
which are not formally defined until the next Chapter. As such, we complete Example 5.4 
at the end of Chapter 6. 



Chapter 6: Affine Approximations 

Remark 6.1: To show that a sub-box of the parameter box W is killed by one of the 
interesting conditions (plus associated killerword) we need to show that at each point in 
the sub-box, the killerword evaluated at that point satisfies the given condition. That is, 
we are simply analyzing a certain function from the sub-box to C. 

As described in Remark 6.5, this analysis can be pulled back from the sub-box in 
question to the unit complex 3-disc A, where A = {(zq, z±, z<l) G C 3 : \zk\ < 1 for k = 
1, 2, 3}. Loosely, we will analyze such a function on A by using Taylor series approximations 
consisting of an affine approximating function together with a bound on the "error" in the 
approximation (this could also be described as a "remainder bound" ) . 
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Problems 6.2: There are two immediate problems likely to arise from this Taylor ap- 
proximation approach. The first problem is the appearance of unpleasant functions such 
as cosh -1 . We have already taken care of this problem by "exponentiating" our prelimi- 
nary parameter space V . This resulted in all functions under consideration being built up 
from the co-ordinate functions V , D' , and R' on V by means of the elementary operations 

Second, for a given "built-up function" the computer needs to be able to compute 
the Taylor approximation, and the error term. This will be handled in this section by 
developing combination formulas for elementary operations (see the Propositions below). 
Specifically, given two Taylor approximations with error terms representing functions g 
and h and an elementary operation on g and h, we will show how to get the Taylor 
approximation with error term for the resultant function from the two original Taylor 
approximations. 

Remark 6.3: We set up the Taylor approximation approach rigorously as follows. The 
notation will be a bit unusual, but we are motivated by a desire to stay close to the notation 
used in the checker computer programs, verify and fudging. However, it should be pointed 
out that the formulas in this Chapter will be superceded by the ones in Chapter 8, which 
incorporate a round-off error analysis. It is the Chapter 8 formulas that are used in verify 
and fudging. 

Definition 6.4: An AffApprox x is a five-tuple (x.f; x.fo, x.fi, x.f 2', x.e), consisting of 
four complex numbers x.f, x.fo, x.fi, x.f 2 and one real number x.e, which represents all 
functions g : A — > C such that 

\g{z Q ,z\, z 2 ) - {x.f + x.fozo + x.fxzx + x.f 2 z 2 )\ < x.e 

for all [zq,z\,z 2 ] E A. That is, x represents all functions from A to C that are x.e— well- 
approximated by the affine function x.f + x.fozo + x.f\Z\ + x.f 2 z 2 . We will denote this 
set of functions associated with x by r(x). 

Remark 6.5: As mentioned in Remark 6.1, given a sub-box to analyze, instead of working 
with functions defined on the sub-box, we will work with corresponding functions defined 
on A. Specifically, rather than build up a function by elementary operations performed on 
the co-ordinate functions L',D',R' restricted to the given sub-box, we will perform the 
elementary operations on the following functions defined on A, 

(Po + ¥3; s + is 3 , 0, 0; 0) (pi + ip 4 ; 0, si + is 4 , 0; 0) (p 2 + ip 5 ; 0, 0, s 2 + is 5 ; 0) 

where (po + ipz,Pi + ip4,p 2 + ip§) is the center of the sub-box in question, and the Sj 
describe the six dimensions of the box. In the computer programs, these three functions 
are called along, ortho, and whirle, respectively, and pi and Sj are denoted pos[i] and size[i], 
respectively. 

After the following Remarks, we state and prove the combination formulas. 

Remarks 6.6: i) We will break with the convention used previously in this paper and 
start the numbering of the Propositions with 6.1. However, we will end this Chapter with 
Example 6.7. 
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ii) The negation of a set of functions is the set consisting of the negatives of the original 
functions, and similarly for other operations. 

Proposition 6.1 (unary minus): If x is the AffApprox (x.f; x.fo, x.fi, x.f2', x.e) then 
— (r(x)) = r(—x) where 

-x = (-x.f; -x.fo, ~x.fi, -x.f 2; x.e). 

Proof: 

\g(z , z x , z 2 ) - (x.f + x.f z + x.fxzx + x.f 2 z 2 )\ < e 

if and only if 

I - g(z , z\, z 2 ) - (-x.f - x.fozo - x.fxzx - x.f 2 z 2 ) \ < e 

I 

Proposition 6.2 (addition): If x is the AffApprox (x.f; x.fo, x.fi, x.f 2', x.e) and y is 
the AffApprox (y.f; y.f , yf u y.f 2 ; y.e), then r(x + y) D r(x) +r(y), where 

x + y = (x.f + y.f; x.fo + y.fo, x.fi + y./i, x.f 2 + y.f 2; x.e + y.e) 

Proof: If g E r(x) and h £ r(y) then we must show that g + h E r(x + y). 

\(g + h)(z ,z u z 2 ) - ((x.f + y.f) + (x.fo + y.fo)zo + (x.fi + y.f\)z x + (x.f 2 + y-f2)z 2 )\ 

< \g(zo, zi, z 2 ) - (x.f + (x.f Q )zo + (x.fi)z-i + (x.f 2 )z 2 ) + 

h(z ,z 1 ,z 2 ) - (y.f + (y.fo)z + (y.fi)zi + (y.f 2 )z 2 )\ 

< \g(zo, z u z 2 ) - (x.f + (x.f )z + (x.fi)z! + (x.f 2 )z 2 )\ + 

\h(z ,z ll z 2 ) - (y.f + (y.f )z + + (y-f2)z 2 )\ 

< x.e + y.e 

I 

We now do subtraction. The statement and proof are essentially the same as for 
addition. The only thing to note is that the errors add. 

Proposition 6.3 (subtraction): If x is the AffApprox (x.f; x.fo, x.fi, x.f 2 ; x.e) and 
y is the AffApprox (y.f; y.f , yfi, y.f 2 ; y.e), then r(x - y) D r(x) - r(y), where 

x-y = (x.f - y.f; x.fo - y.fo, x.f\ - y.fx, x.f 2 - y.f 2 ; x.e + y.e) 
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We now state variations on Propositions 6.2 and 6.3 whose usefulness will not be 
apparent until Chapter 8, when we incorporate round-off error into these formulas. In 
what follows, a "double" corresponds to a real number, and has an associated AffApprox, 
with last four entries zero. 

Proposition 6.4 (addition of an AffApprox and a double): If x is the AffApprox 
(x.f; x.fo, x.fi, x.f 2 ; x.e) and y is a double, then r(x + y) D r(x) +r(y), where 

x + y = (x.f + y; x.fo, x.fi, x.f 2 ; x.e). 



Proposition 6.5 (subtraction of a double from an AffApprox): If x is the AffApprox 
(x.f; x.fo, x.fi, x.fi', x.e) and y is a double, then r(x — y) D r(x) — r(y), where 

x-y= (x.f - y; x.fo, x.fi, x.f 2 ; x.e). 



Proposition 6.6 (multiplication): If x is the AffApprox (x.f; x.fo, x.fi, x.f 2 ; x.e) 
and y is the AffApprox (y.f; y.f , yfi, y.f 2 ; y.e), then r(x x y) D r(x) x r(y), where 

xxy = (x.f x y.f; x.f x y.f + x.fo x y.f, x.f x y.fi + x.fi x y.f, x.f x y.f 2 + x.f 2 x y.f; 

(size(x) + x.e) x (size(y) + y.e) + (\x.f\ x y.e + x.e x \y.f\)) 

with size(x) = \x.f \ + \x.fi\ + \x.f 2 \ and size(y) = \y.f \ + \y.fi\ + \y.f 2 \ 

Proof: If g G r(x) and h G r(y) then we must show that g x h G r(x x y). That is, we 
need to show 

\(g x h)(zo,Zi,z 2 ) - ((x.f x y.f) + 
(x.f x y.fo+x.fo x y.f)z + (x.f x y.fi + x.fi x y.f)z x + (x.f x y.f 2 + x.f 2 x y.f)z 2 )\ 

< (size(x) + x.e) x (size(y) + y.e) + (\x.f\ x y.e + x.e x \y.f\) 

Note that for any point (zq,Zi,z 2 ) G A and any functions g G r(x) and h G r(y) we can 
find complex numbers u, v with \u\ < 1 and \v\ < 1, such that 

g(zo, zi,z 2 ) = x.f + (x.fozo + x.f\z x + x.f 2 z 2 ) + (x.e)u 

and 

h(zo, zi, z 2 ) = y.f + (y.fozo + y.fizi + y.f 2 z 2 ) + (y.e)v. 
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Multiplying out, we see that 

(g x h)(z , z x , z 2 ) = (x.f x y.f)+ 
(x.f x y.f + x.fo x y.f)z + (x.f x y.f x + x.f x x y.f)z x + (x.f x y./ 2 + x.f 2 x y.f)z 2 + 

(x.f x y.e)v + (x.e x y.f)u+ 
((x.f z + x.f x z x + x.f 2 z 2 ) + (x.e)u) x ((y.f z + y.fiZi + y.f 2 z 2 ) + (y.e)v) 

Hence, 

\(g x h)(z ,z x ,z 2 ) - ((x.f x y.f) + 
((x.f xy.f + x.f x y.f)z + (x.f x y.f x + x.f x x y.f)z x + (x.f x y.f 2 + x.f 2 x y.f)z 2 ))\ 

< (\x.f\y.e + x.e\y.f\) + (size(x) + x.e) x (size(y) + y.e). 



Proposition 6.7 (an AffApprox multiplied by a double): If x is the AffApprox 
(x.f; x.fo, x.f x , x.f 2 ; x.e) and y is a double, then r(x xi/)D r(x) x r(y), where 

xxy = (x.f x y; x.fo x y, x./i x y, x./ 2 x x.e x |y|) 



Proposition 6.8 (division): If x is the AffApprox (x.f; x.fo, x.f x , x.f 2 ; x.e) and y is 
the AffApprox (y.f; y.f , yf x , y.f 2 ; y.e), then r(x/y) D r(x)/r(y), where 

x/y = (x.f/y.f; (-x.f x y.f + x.fo x y.f)/(y.f 2 ), (-x.f x y./x + x./i x y.f)/(y.f), 

(-x.f x y./ 2 + x./ 2 x y.f)/(y.f); 
(\x.f\ + size(x) + x.e)/(\y.f\ - (size(y) + y.e))- 

((\x.f\/\y.f\ + stze(x)/\y.f\) + \x.f\size(y)/(\y.f\\y.f\))) 

Of course, we have to avoid division by zero. That is, we demand that \y.f \ > size(y)+y.e. 

Proof: For notational convenience, denote (x.foZo + x.f x z x + x.f 2 z 2 ) by x.fkZk and sim- 
ilarly for y.fkZk and so on. As above, note that for any point (zo, z x , z 2 ) G A and any 
functions g G r(x) and h G r(y) we can find complex numbers u, v with \u\ < 1 and 
\v\ < 1, such that 

g(z , z x , z 2 ) = x.f + (x.f k z k ) + (x.e)u 

and 

h(z , zx, z 2 ) = y.f + (y.fkZk) + (y.e)v. 

We compare (g/h)(zo, z x ,z 2 ) with its putative affine approximation. That is, we ana- 
lyze 

\(x.f + (x.fkZk) + (x.e)u)/(y.f + (y.f k z k ) + (y-e)v)- 

y * j 
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Putting this over a common denominator of \(y.f 2 )(y.f + (y.fkZk) + (y.e)v)\ and cancelling 
equal terms (in the numerator) we are left with a quotient whose numerator is 

\x.e(y.f 2 )u - (x.f k )y.f(y.f k )z k - x.f(y.f k 2 )z k + 

(x-f)y.f(y.e)v + x.fk(y-f)y-e(v)z k - x.f(y.f k )y.e(v)z k \. 

We must show this (first) quotient is bounded by 

(\x.f\ + size(x) + x.e)/(\y.f\ - (size(y) + y.e))- 

((\x.f\/\y.f\ + size(x)/\y.f\) + \x.f\size(y)/(\y.f\\y.f\)). 

Putting this over a common denominator of \y-f\ 2 (\y-f\ — (size(y) + y.e)) and cancelling 
equal terms (in the numerator) we are left with a second quotient, whose numerator is 

x.e\y.f\ 2 - (-\x.f\\y.f\y.e - size(x)\y.f\(size(y) + y.e) - \x.f\size(y)(size(y) + y.e)) 

and we see that all terms in this numerator are positive. Further, the terms in the nu- 
merators of the first and second quotients correspond in a natural way, and each term in 
the numerator of the second quotient is greater than or equal to the absolute value of its 
corresponding term in the numerator of the first quotient. 

Finally, because the denominator in the second quotient is less than or equal to the 
absolute value of the denominator in the first quotient, we see that the absolute value of 
the first quotient is less than or equal to the second quotient, as desired. I 

We present a couple of variations on Proposition 6.8 which will be useful when we do 
round-off error. 

Proposition 6.9 (division of a double by an AffApprox): If x is a double and y is 

the AffApprox (y.f; y.f , yf u y.f 2 ; y.e), then r(x/y) D r(x)/r(y), where 

x/y = (x/y.f; -x x y.f /(y.f 2 ), -x.f x y.fi/(y.f), -x.f x y.f 2 j(y.f\ 

{\Ali\v-f \ ~ (size(y) +y.e))- 

(\x\/\y.f\ + \x\size(y)/(\y.f\\y.f\))) 

Of course, we have to avoid division by zero. That is, we demand that \y.f\ > size(y)+y.e. 
I 

Proposition 6.10 (division of an AffApprox by a double): If x is the AffApprox 
(x.f; x.fo, x.f x.f 2, x.e) and y is a double, then r(x/y) D r(x)/r(y), where 

x/y = (x.f/y; x.f /y, x.J \/y, x.f 2 /y; x.e/\y\) 

Of course, we have to avoid division by zero. That is, we demand that \y\ > 0. I 

Finally, we do the square root. 
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Proposition 6.11 (square root): If x is the Aff Approx (x.f; x.fo, x.fi, x.fi\ x.e) then 
r(\fx) D yjr{x), where 

r ( I — 1 x -fo x -h x -f2 



2v^7' 2v^7' 2v^7' 

VWT\ - + y/fcf\ - (size(x) + x.e))) 

2VF-/I 

Of course, we have to avoid division by zero, and taking the real square root of a negative 
number. In particular, when \x.f\ > size(x) + x.e we get the formula above, but if this 
does not hold, then we use the crude estimate (0; 0, 0, 0; ^/\x.f\ + size(x) + x.e). 

Proof: As above, note that for any point (20,-21,-22) £ A and any function g G r(x) we 
can find a complex number u with \u\ < 1, such that 

g(z , zi,z 2 ) = x.f + (x.f k z k ) + (x.e)u. 

Also, because \x.f\ > size(x) + x.e. we see that the argument of x.f + (x.fkZk) + (x.e)u is 
within 7r/2 of the argument of x.f, and therefore, we can require that >Jg(zo, z±, Z2) has 
argument within w/4 of the argument of \Jx.f . 
We need to show that 

\y/x.f + x.f k z k + (x.e)u - (a/x7/+ X 'JjL 7 )z k \ 



< V\ x -f\ - ( ~ /;— T7 + V\ x -f\ - (size(x) +x.e) ) 
2VF-/I 

Or, after multiplying both sides by yj\x.f\, 

V x.f (x.f + x.fkZk + (x.e)u) - (x.f + (x.f k )z k /2) 



< (\x.f\ - size(x)/2) - ^\x.f\(\x.f\-(size(x) + x.e)) 

The two sides of the inequality are of the form A — B and C — D, and we "simplify" by 
multiplying by -jfqr§ and We now show that the (absolute value of the) left-hand 

numerator is less than or equal to the right-hand numerator. Later, we will show that 
the (absolute value of the) left-hand denominator is larger than or equal to the right-hand 
denominator. The left-hand numerator is 

\x.f(x.f + x.fkZk + (x.e)u) - (x.f + (x.f k )z k /2) 2 \ 

= \x.f 2 + x.f(x.fk)z k + x.f(x.e)u - x.f 2 - x.f(x.f k )z k - (x.fl)z 2 k /A\ 

= \x.f(x.e)u-(x.f 2 )z 2 /A\ 
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The right-hand numerator is 

(\x.f \ — size(x) /2) 2 — — (size(x) + x.e)) 

= \x.f\ 2 — \x.f\size(x) + size(x) 2 /4 — \x.f\ 2 + \x.f\size(x) + \x.f\x.e 

= \x.f\x.e + size(x) 2 /4 

So the left-hand numerator is indeed less than or equal to the right-hand numerator. 

We now compare the denominators, but only after dividing each by y/\x.f\. The 
left-hand denominator is 

\y/x.f + x.f k z k + (x.e)u+ (\/xJ + ' 



2v^7 

while the right-hand denominator is 

V\x-f\ ~ ^/r-rr + V\x-f\ - (size(x) + x.e) 
2VF-/I 

The claim that the left-hand denominator is greater than or equal to the right-hand de- 
nominator is a bit complicated. First, compare the y/x.f term and the y/\x.f\ terms. 
They are the same distance from the origin. Next, note that as z& and u take on all values, 
x.f + x.fkZk + (x.e)u describes a disk centered at x.f and whose radius is less than \y/x.f\. 
Hence, \Jx.f + x.fkZk + (x.e)u describes a convex symmetric (about the line formed by 
the origin and x.f) set centered at y/x.f. Further, y/x.f + y/x.f + x.fkZk + (x.e)u de- 
scribes a convex symmetric (about the line formed by the origin and x.f) set centered at 
2y/x.f. In any case, it is easy enough to see that no points on this convex symmetric set 
get closer to the origin than y/\x.f\ + y/\x.f\ — (size(x) +x.e). 
Finally, because \ x ^pL \ < siz < x ) no po ints of 

y/xj+ y/ x.f + x.f k z k + (x.e)u + X J^X 

2y/x.f 

can get closer to the origin than 

V\xT\ + V\ x -f\ - (size(x) + x.e) - ^== 



Example 6.7 (Continuation of Example 5.4): We repeat the description of the sub- 
box under investigation. 

The sub-box Z(s01011) with 

s = 001000110001110111001111000101111111101111100111001111000001111011110111 
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is the region where 



f -1.381589027741 . 
-1.378124546093. 

0.999893182771 . 
-2.535837191243. 
2.535404997792 . . 
V -0.001953125000 



. . < Re(L') < -1.379848991182 . 
. . < Re(D') < -1.376574349753. 
. . < Re(R') < 1.001274250703 . . 
. . < Im(L') < -2.534606799593 . 
. < Im(D') < -2.534308843448 . 
. . . < Im(R') < 0.000000000000 . 



For this sub-box, we get 

( -0.8677851121 + il. 4607429651; ^ 
0.0000248810 - z0.0003125810, 
0.0000000000 + zO.0000000000, 
0.0000000000 + zO.0000000000; 
V 0.0000000289 / 

( 0.0000000000 + zO.0000000000; ^ 
0.0000000000 + zO.0000000000, 
0.0000000000 + zO.0000000000, 
0.0000000000 + zO.0000000000; 
V 0.0000000000 / 



/ = 



and 



w = 



( -0.5845111829 + z0.4773282853; \_ 
0.0000000000 + zO.0000000000, 
-0.0000296707 - z0.0001657332, 
-0.0004345111 - zO.0001209539; 
0.0000002590 
( -0.2832291572 + z0.9833572297; \ 
0.0000000000 + zO.0000000000, 
0.0000515806 - z0.0001129408, 
-0.0005778031 + z0.0002005440; 
0.0000002806 



/ 0.0000000000 + zO.0000000000; \ 
0.0000000000 + zO.0000000000, 
0.0000000000 + zO.0000000000, 
0.0000000000 + zO.0000000000; 
V 0.0000000000 J 

( -0.3006023265 - z0.5060039953; \ 
' -0.0000909686 - z0.0000593570, 
0.0000000000 + zO.0000000000, 
0.0000000000 + zO.0000000000; 
V 0.0000000301 / 



( -0.2840228472 + z0.9825063583; \ 
0.0000000000 + zO.0000000000, 
0.0000516606 - z0.0001128245, 
0.0005776611 - z0.0001998632; 
0.0000006462 
( -0.5846352333 + z0.4764792236; \ 
0.0000000000 + zO.0000000000, 
-0.0000294917 - z0.0001656653 
0.0004341392 + i0.0001213070; 
V 0.0000005286 



calculating g = f 1 wf 1 w 1 / w 1 fw 1 / w 1 / wf 1 wfww gives 



9 = 



( -0.5764337542 + i0.4752708071; ^ 
-0.0031657223 - z0.0001436786 
-0.0017723577 + z0.0000352928 
-0.0011623491 + z0.0017516088 

V 0.0008229225 / 
f -0.2861207992 + z0.9766064999; \_ 

-0.0002777968 + i0.0020330488, 
0.0000837571 + z0.0010241875, 
0.0028322367 - i0.0005972336; 

V 0.0018172437 / 



( -0.2704033973 + z0.9822741250; ^ 
-0.0045902952 - z0.0019135041, 
-0.0026219461 - z0.0007506230, 
-0.0002823450 + z0.0033805602; 

V 0.0008037640 / 
( -0.5861133046 + z0.4624368851; \_ 

-0.0021932627 + i0.0040523411 : 
-0.0008612361 + z0.0022394639 
0.0061581377 - i0.0005862070; 

V 0.0017738513 / 
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We then get 



length(g) 



( -1.3588762105 - z2.4897230182; \_ 
0.0030210500 - z0.0182284729, 
0.0007938572 - z0.0096614614, 
-0.0122034521 + z0.0074353043; 
V 0.0080071969 / 



and 



length(g) 



( 0.9825397896 - z0.0008933519; ^ 
0.0053701602 + z0.0037789019, 
0.0028076072 + iO.0018421952, 
-0.0002400615 - i0.0049443045; 
V 0.0027802966 / 



This isn't quite good enough to kill the sub-box, since \length(g)/L'\ can be high as 
1.0001951323. 

When we subdivide Z(s01011), we have to analyze two sub-boxes, Z(sOlOllO) and 
Z(sOlOlll). For Z(sOlOllO), the same calculation on the region 
-1.381589027741073400 < Re(L') < -1.379848991182205200 
-1.378124546093485700 < Re(D') < -1.376574349753672900 
0.999893182771602220 < Re(R') < 1.001274250703607400 
-2.535837191243490300 < Im(L') < -2.534606799593201600 
-2.535404997792558600 < Im(D') < -2.534308843448505900 
-0.001953125000000000 < Im(R') < -0.000976562500000000 
gives 

/ 0.9814518667 + z0.0008103446; \ 
0.0053616729 + i0.0037834001, 
0.0028027236 + z0.0018435245, 
-0.0013175066 - z0.0032448794; 
V 0.0019033926 / 



length(g) 



and we can then bound | len 9^{g) | < 0.9967745579, which kills Z(sOlOllO). 
On Z(sOlOlll), the calculation gives 



length(g) 



/ 0.9836225919 - z0.0025990177; \ 
0.0053786346 + i0.0037743930, 
0.0028124892 + z0.0018408583, 
-0.0013333182 - i0.0032343347; 

V 0.0019044429 / 



and | iew y g> | < 0.9989610507, which kills Z(sOlOlll). 
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Chapter 7: Complex Numbers with Round-Off Error 

Remark 7.1: The theoretical method for proving Theorem 0.2 has been implemented 
on the computer programs verify and fudging given in the Appendices. To make this 
computer-aided proof rigorous, we need to deal with round-off error in calculations. 

One approach to round-off error would be to use interval arithmetic packages to carry 
out all calculations with floating-point numbers (also called "doubles"), or to generate our 
own version of these packages. However, it appears that this would be much too slow given 
the size of our collection of sub-boxes and conditions and killer words. 

To solve this problem of speed, we implement round-off error at a higher level of 
programing. That is, we incorporate round-off error directly into AffApprox's. This ne- 
cessitates that we incorporate it into complex numbers as well. 

Definition 7.2: In the next Chapter we work with AffApprox's. In this section we show 
how to do standard operations on complex numbers while keeping track of round-off error. 
There are two types of complex numbers to consider: 

1. ) An XComplex corresponds to a complex number that is represented exactly. Thus, 
it simply consists of a real part and an imaginary part. 

2. ) An AComplex corresponds to an "interval" that contains the complex number in 
question. Thus, it consists of an XComplex and a floating-point number representing the 
error. In particular, (x;e) represents the set of complex numbers {w : \w — x\ < e }. 
Following the notation of Chapter 6, we could now denote this set of complex numbers 
by r(x, e), but instead we suppress mention of the r functions throughout the rest of this 
section. It seems preferable to abuse notation in this fashion in the interests of simplicity. 

Remark 7.3: In general, our operations act on XComplexes and produce AComplexes, 
or they act on AComplexes and produce AComplexes . In one case, the unary minus, an 
XComplex goes to an XComplex. In the calculations that follow the effect on the error is 
the whole point. 

Conventions 7.4: We begin, by writing down our basic rules, which follow easily from 
the IEEE-754 double-precision standard for machine arithmetic. (Actually, the "hypot" 
function h(a, b), which computes by elaborate chicanery Va 2 + b 2 , is not part of the IEEE 
standard, but satisfies the appropriate standard according to the documentation provided 
(see [Kahan]).) The operations here are on double-precision floating-point real numbers 
(doubles) and we denote a true operation by the usual symbol and the associated machine 
operation by the same symbol in a circle, with two exceptions: a machine square root ^/a is 
denoted \fa and the machine version of the hypot function is denoted h . Perhaps a third 
exception is our occasional notation of true multiplication by the absence of a symbol. 
The standard number used in error analysis is an "EPS." It depends on the number of bits 
used to store floating-point numbers, and this can vary from machine to machine. Finally, 
it should be noted that our analysis breaks down when there is "underflow," so in the 
computer programs we ask the computer to inform us if underflow has occurred. 
As in Chapter 6, we now break with the usual numbering convention. 

Basic Properties 7.0 (assuming no underflow): 
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1 + k x EPS = 1 © (k <8 EPS) when is an integer which is not huge in absolute 
value. 

2 fc x A = 2 k © A when /c is an integer, 2 k A is not infinity. 



(a 


+ b) - 

T V) 


- (a 


ffi Ml < (EPS/2)\a 


4- 61 

~ 1 


(a 


+ &)" 


- (a 


mM|< f£'P t S72 N |l<7 


ffi 61 


(a 


-6)- 


- (a 


qM < (EPS/2)\a 


- 61 


(a 


-6)- 


- (a 


0o) < (EPS/2) a 


^ L 1 

e 6| 


(a 


x 6) - 


-(a 


(8)6)| < {EPS/2)\a 


x 6| 


(a 


x 6) - 


- (a 


(8)6)| < (EPS/2)\a 


©6| 


l( 


a/6) - 


- (a 


6)| < (EPS/2)\a 


/6| 


|(a/6) - 


(aC 


3 6)| < (EPS/2)\a(. 


2>6| 




\\fa 


V 


^| < (£PS/2)|<y/a| 






\y/a 


- V 


^| < (^5/2)1^1 





|/i(o,6) - /i (a,6)| < (EPS)\h(a,b)\ 
\h(a,b) - h (a,b)\ < (EPS)\h (a,b)\ 

From these formulas, we immediately compute the following. 



(1- 


EPS/2)\a + b\ < 


a 


© 6| 


< 


(1 + EPS/2) 


a + b\ 


(1- 


EPS/2)\a®b\ < 


a 


+ 6| 


< 


(1 + EPS/2) 


a©6| 


(1- 


EPS/2)\a-b\ < 


a 


e6| 


< 


(1 + EPS/2) 


a-b\ 


(1- 


EPS/2)\aQb\ < 


a 


-6| 


< 


(1 + EPS/2) 


a©6| 


(1- 


EPS/2)\a x 6| < 


a 


©6| 


< 


(1 + EPS/2) 


a x 6 


(1- 


EPS/2)\a(g)b\ < 


a 


x 6| 


< 


(1 + EPS/2) 


a © 6 


(1 


-EPS/2)\a/b\ < 


a 


6| 


< 


(1 + EPS/2) 


a/6| 



(1- EPS/2)\a0b\ < \a/b\ < (1 + EPS/2)\a 6| 
(1 - £PS/2) | Va| < | < (1 + EPS/2) | 
(1-£?P5/2)|^| < \y/E\ < (l + EPS/2)\^\ 
(1- EPS)\h(a,b)\ < \h (a,b)\ < (1 + EPS)\h(a,b)\ 
(1- EPS)\h (a,b)\ < \h(a,b)\ < (1 + EPS)\h (a, b)\ 
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Of course, we can also get the following type of formula, which is sometimes convenient, 
for example, in the proof of Lemma 7.2. 

( - , 1 -EPs )\a®b\ < \a + b\ < ( - 1 - MEl )\a®b\ 

1 + 2 1 2 

Proof: There is a finite set of numbers which are representable on the computer, and the 
IEEE standard states that the result of an operation is always the closest representable 
number to the true solution. Ignoring technicalities, a non-zero floating-point number is 
represented by a fixed number of bits of which the first determines the sign of the number, 
the next m represent the exponent, and the remaining n represent the mantissa of the 
number. Because our non-zero numbers start with a 1, that means the n mantissa bits 
actually represent the next n binary digits after the 1. That is, the mantissa is actually 
1.&1&2&3--&71- With this divyying up of the m + n + 1 bits among exponent and mantissa, 
EPS would be 2~ n and EPS/2 would be 2~ ( - n + 1 \ 
Given this set-up, properties of the form 

\(a + b) - (a ©6)| < (EPS/2)\a + b\ 

follow immediately. Then, 

\(a + b) - (a ©6)| < (EPS/2)\a®b\ 

follows because the true answer has "exponent" which is less than or equal to the exponent 
of the machine answer. | 

Before starting in with our Propositions, we prove a couple of lemmas. 

Lemma 7.0: 

(1 - EPS) © \a © b\ < \a + b\ < (1 + EPS) © \a © b\ 
Analagous formulas hold for — , *, /, ^J~. 

Proof: Assume a + b > 0. If (1 + EPS) © (a © b) < (a + b) then the machine 
number (1 + EPS) © (a © b) is a better approximation to a + b than a © 6, because 
(a © b) < (1 + EPS) © (a © b). This contradicts the IEEE standard. The case a + b < 
can be handled similarly, and the case a + b = is trivial. Similarly for the left-hand 
inequality. I 

Lemma 7.1: 

(1 + EPS/2) a A < (1 + kEPS) © A 

where A > 0, and a and k are (not huge) integers, such that for a even, k = | + 1 and for 
a odd, k = s±i + 1. 
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Proof: 

(1 + EPS/2) a A < (1 - EPS/2)(\ + kEPS)A < (1 + kEPS) <g> A 

The first inequality holds if o and are as in the Lemma, and the second inequality is a 
consequence of one of the formulas preceding Lemma 7.0 (A > 0). I 

We start the operations. We will give proofs for most, the others should be straight- 
forward to derive. Note that for an XComplex, x = (x.re,x.im), and for an AComplex, 
x = (x.re,x.im,x.e). 

Proposition 7.1 (-X): 

— x = (— x.re, — x.im). | 

Proposition 7.2 (X + D) (an XComplex and a double added together, which yields an 
AComplex) : 

x + d= (x.re © d, x.im; (EPS/2) © \x.re © d\) 

Proof: The error is bounded by 

| (x.re + d) — (x.re © d)\ 

< (EPS/2)\x.re®d\ 
= (EPS/2) © \x.re®d\ 

I 

Proposition 7.3 (X - D) (a double subtracted from an XComplex, which yields an 
AComplex) : 

x - d = (x.re © d, x.im; (EPS/2) © \x.re © d\) | 

Proposition 7.4 (X + X) (an XComplex and an XComplex added together, which yields 
an AComplex): 

x+y = (x.reQ)y.re,x.im(By.im; (EPS /2)®((l J r EPS)®(\x.re@y.re\@\x.im®y.im\))) 

Proof: The error is bounded by 

| (x.re + y.re) — (x.re © y.re) \ + \ (x.im + y.im) — (x.im © y.im) \ 

< (EPS/2)(\x.re © y.re\ + \x.im © y.im\) 
< (EPS/2)((1 + EPS) © (\x.re © y.re\ © \x.im © y.im\)) 
= (EPS/2) © ((1 + EPS) © (\x.re © y.re\ © \x.im © y.im\)) 
To go from line 2 to line 3 we used Lemma 7.0. I 
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Proposition 7.5 (X - X): 

x+y = (x.reQy.re,x.imQy.im; (EPS /2)®((l J r EPS)®(\x.reQy.re\@\x.imQy.im\))) 

I 

Proposition 7.6 (A + A) (an AComplex and an AComplex added together, which yields 
an AComplex): 

x + y = (re, im; e) where 

re = x.re © y.re 

im = x.im © y.im 

e = (1 + 2EPS) © (((EPS/2) © (\re\ © \im\)) © (x.e © y.e)) 

Proof: The error is bounded by the sum of the contributions from the real part, the 
imaginary part, and the two individual errors: 

| (x.re © y.re) — (x.re + y.re)\ + \(x.im © y.im) — (x.im + y.im)\ + (x.e + y.e). 

< (EPS/2)\x.re © y.re\ + (EPS/2)\x.im © y.im\ + (1 + EPS/2)(x.e © y.e) 

< (1 + EPS/2)(EPS/2)(\x.re © y.re\ © \x.im © y.im\) + (1 + EPS/2)(x.e © y.e) 
= (1 + J EP5 , /2)(( J EP5'/2)(|a:.re © y.re| © |x.im © y.im\) + (x.e © y.e)) 
< (1 + EPS/2) 2 (((EPS/2)(\x.re © y.re| © \x.im © y.im|)) © (x.e © y.e)) 

< (1 + 2EPS) © (((EPS/2) © flx.re © y.re| © |x.im © y.im|)) © (x.e © y.e)) 

I 

The hierarchy for machine operations is the same as that for true operations, so one 
pair of parentheses is unnecessary and will often be omitted in what follows. 

Proposition 7.7 (A - A): 

x + y = (re, im; e) where 
re = x.re y.re 
im = x.im y.im 

e = (1 + 2EPS) © (((EPS/2) © (|re| © \im\)) © (x.e © y.e)) I 

Proposition 7.8 (X x D): 

x x d = (re, im; e) where 
re = x.re © y 
im = x.im © y 

e = (EPS/2) © ((1 + EPS) © (|re| © \im\)) 
Proof: The error is bounded by 

| (x.re x y) — (x.re <g> y)\ + | (x.im x y) — (x.im © y) | 
< (EP5/2)|x.re©y| + (£PS/2)|a:im©y| 
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= (EPS/2)(\x.re y\ + \x.im y|) 
< {EPS/2) ((1 + EPS) (\x.re y\ © \x.im y|)) 

I 

Proposition 7.9 (X / D): 

x/d= (re, im; e)where 
re = x.re y 
im = x.im y 

e = (EPS/2) ((1 + EPS) (|re| © |im|)) I 

Proposition 7.10 (X x X): 

x x y = (re, im; e) where 

re = rel re2, with rel = x.re y.re and re2 = x.im y.im 
im = iml © im2, with iml = x.re y.im and im2 = x.im y.re 
e = EPS ((1 + 2EPS) ((|rel| © |re2|) © (\iml\ © \im2\))) 

Proof: The error is bounded by the sum of the contributions from the real part and 
the imaginary part: 

| (x.re x y.re — x.im x y.im) — ((x.re y.re) (x.im y.im))\ 

+ \(x.re x y.im + x.im x y.re) — ((x.re y.re) © (x.im y.im))\ 

We want to bound this by a machine formula. Let's begin by bounding 

\(x.re x y.re — x.im x y.im) — ((x.re y.re) (x.im y.zm))| 

by a machine formula. 

| (x.re x y.re — x.im x y.zm) — ((x.re y.re) (x.im y.zm))| 

< | ((x.re x y.re) — (x.im x y.im)) — ((x.re y.re) — (x.im y.zm))| 

+ | ((x.re y.re) — (x.im y.im)) — ((x.re y.re) (x.im y.zm))| 

< | (x.re x y.re) — (x.re y.re)| + |(x.zm x y.im) — (x.im y.im)\ 

+ (EPS/2)\(x.re y.re) - (x.im yim)| 

< (EPS/2) | (x.re y.re) | + (EPS/2) | (x.im y.im) | 

+ (EPS/2)(\x.re y.re| + |x.im y.im|) 

= (EPS/2)(2)(\x.re y.re| + |x.im y.im\) 
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< EPS(1 + EPS/2)(\x.re © y.re\ © \(x.im © y.im\) 

Almost the exact same calculation produces the analagous formula for the imaginary 
contribution, and we now combine the two to get a bound on the total error. 



< EPS(l + EPS/2)(\x.re © y.re\ © \x.im © y.im\) 

+ EPS(1 + EPS/2)(\x.re © y.im\ © \x.im © y.re\) 

= EPS(1 + EPS/2)((x.re © y.re\ © \x.im © y.im\) 

+ (\x.re © y.im\ © \x.im © y.re\)) 

< EPS(1 + EPS/2) 2 ((\x. re © y.re\ © \x.im © y.im\) 

© (|x.re © y.im\ © |xim © y.re\)) 

< EPS © ((1 + 2EPS) © (flx.re © y.re\ © |x.im © y.im|) 

© (\x.re © j/.zm| © \x.im © j/.re|))) 



Proposition 7.11 (D /X): 

x/y = (re, im; e) where 

re = (x © j/. re) nrm where nrm = y.re © y.re © yim © y.im 

im = — (x © y.im) nrm 

e = (2EPS) © ((1 + 2EPS) © (|re| © |zm|)) 

Proof: The true version of x/y is equal to (x x y.re + i(— x x y.im)) /(y.re 2 +y.im 2 ) 
and we need to compare this with the machine version to find the error. Further, this 
error is less than or equal to the sum of the real error and the imaginary error. Thus, we 
start with the real calculation (as above, we use nrm to represent the machine version of 
y.re 2 + y.im 2 ). 

x x y.re 

2 : — 2 — {{x © y.re) nrm)\ 



y.re z + y.im z 

x®y.re. .x®y.re xxy.re. .x x y.re x x y.re 
< \(x © y.re) nrm 1 + | 1 + 



nrm 1 1 nrm nrm 1 1 nrm y.re 2 + y.im 2 
Before continuing, let's compare and tt — : — 2 by developing a formula for 

03 i nrm y .re z -\-y .im^ J e 



comparing a2 ^_ b2 and its associated 
Lemma 7.2: 



< (EPS + (EPS/2) 2 )- J 



nrm a 2 + 6 2 //r//; 
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where nrm = a ® a ®b ®b. 
Proof: We compute that 



) 2 x nrm < a 2 + b 2 < (- J" ) 2 x nrm, 



v l + EPS/2' ~ - v l- EPS/2 

hence 



^— (1 - EPS/2) 2 < . 1 I0 < (1 + EPS/2) 2 , 
nrm a 2 + nrm 



It then follows that 



' 1 1 < —(1 + EPS/2) 2 ' 



nrm a 2 + b 2 nrm nrm 
'" -((l + £PS/2) 2 -l) = (£PS+(£PS/2) 2 )- J: 



nrm nrm 



Getting back to our main calculation (with nrm = y.re y.re © yim y.im), 

.. . x®y.re. .x®y.re xxy.re. .xxy.re x x y.re 

| (x y.re) nrm 1 + | 1 + | 5 : — ~ 

nrm nrm nrm nrm y.re z + y.im A 

< (EPS/2) lx ® y - rel + (EPS/2) lx ® y - rel + (EPS+(EPS/2) 2 ) lxXy - rel 
nrm nrm nrm 

= (EPS/2)(-^—)(2\x y.rel + (2 + EPS/2) x \x x y.re\) 
nrm 

< (EPS/2)(-^—)(2\x y.re| + (2 + EPS/2)(1 + EPS/2) x \x y.re\) 

= (EPS/2)(—)(\x y.re\){2 + (2 + EPS/2)(1 + EPS/2)) 
nrm 

< (EPS/2)(4 + 3EPS/2 + (EPS/2) 2 )(\x y. re |)(— ) 

nrm 

< (EPS/2)(4 + 3EPS/2 + (EPS/2) 2 )(1 + EPS/2)(\x y.re\ nrm) 

< {2EPS){\ + 3EPS/8 + (EPS/4) 2 )(1 + EPS/2)(\(x y.re nrm)\) 

We also get the analagous formula for the imaginary contribution for the error, so our total 
error is bounded by 

(2EPS)(l + 3EPS/8 + (EPS/4) 2 )(l+EPS/2)((\(x®y.re)0nrm\) + (\(x®y.im)0nrm\)) 

< (2EPS)(l+3EPS/8+(EPS/4) 2 )(l+EPS/2) 2 ((\(x®y.re)0nrm\)®(\(x®y.im)0nrm\)) 
< (2EPS)(1 - EPS '/2)(1 + 2EPS)((\{x y.re) nrm|) © (\(x yim) nrm|)) 
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< (2EPS) © ((1 + 2EPS) © {{\{x®y.re) 0nrm\) © (\(x®y.im) nrm\))) 
Here we used the fact that 

(1 + 3EPS/8 + (EPS '/4) 2 )(1 + EPS/2) 2 < (1 - EPS/2) (1 + 2EPS) 



Proposition 7.12 (X / X): 

x/y = (re,im;e) where 

re = (x.re © y.re © x.im © y.im) nrm where nrm = y.re © y.re © yim © y.im 
im = (x.im © y.re x.re © y.zm) nrm 
e = (5PPS/2) © ((1 + 3PPS) © A) where 

A = ((|x.re © y.re\ © jx.zm © y.im\) © (|xim © y.re| © \x.re © y.im|)) nrm 
Proof: 

It will be useful to begin by comparing (a x b + c x d) / (e x e + f x f) with (a © 6 © 
c © cf) nrm. 

Lemma 7.3: 

| (a © b © c © d) nrm - (axf)+cxd)/exe + /x/)| 

< (EPS/2)(\ + EPS/2)(5 + 3EPS/2 + (EPS/2) 2 )(\a © 6| © |c © d\)/nrm 

Proof: At one point in the proof, we will use the fact that |A©P|<|A|©|P|. This 
fact is easily seen by considering the various possibilities for the sign of A versus the sign 
of B. At a different point in the proof, Lemma 7.2 with a, b replaced by e, / will be used. 

| (a © b © c © d) nrm - (a x b + c x d)/e x e + / x /) | 

< |(a©6©c©(i)0 nrm — (a © b © c © d)/nrm\ 

+ | (a © b © c © d) /nrm -(ax6 + cxd)/(exe + /x/)| 

< (EPS/2) | (a © 6 © c © <i)/nrm| + |(a © b © c © d)/nrm — (a © 6 + c © <i)/nrm| 

+ | (a © 6 + c © d) /nrm — (a x b + c x d) / (e x e + f x f)\ 



< (EPS/2)\a®b® c© d\/nrm+ \ (a <g> b ® c <g> d) - (a ©6 + c© d)\/nrm 
+ | (a © b + c © d)/nrm — (a x b + c x d)/nrm\ 

+ \ (a x b + c x d) /nrm - (a x b + c x d)/(e x e + f x f)\ 
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< (EPS/2)\a®b® c® d\/nrm + (EPS/2)\a®b ® c® d\/nrm 

+ | (a ® b + c ® d)/nrm — (a x b + c x d)/nrm\ 

+ \ (a x b + c x d) /nrm - (a x b + c x d)/(e x e + f x f)\ 



< (EPS/2)\a®b@ c® d\/nrm + (EPS/2)\a <g> b © c® d\/nrm 

+ | (a ® b — a x b)/nrm — (c® d — c x d)/nrm\ 

+ (EPS + (EPS/2) 2 )\(a xb + cx d)\/nrm 

< (EPS)\a ® b © c ® d\/nrm + (EPS/2)\a ® b\/nrm + (EPS/2)\c® d\/nrm 

+ (EPS + (EPS /2) 2 ){l + EPS/2)\(a ®b + c® d)\/nrm 

< (EPS)\a®b®c®d\/nrm+ (EPS/2)\a ® b\/nrm + (EPS/2)\c® d\/nrm 

+ (EPS + (EPS '/2) 2 )(1 + EPS/2) (| a ®b\ + \c® d\)/nrm 

< (EPS)\a®b® c® d\/nrm 

+ (EPS/2 + (EPS + (EPS/2) 2 )(l + EPS/2))(\a ®b\ + \c® d\)/nrm 

< (EPS)\a®b® c® d\/nrm 

+ (EPS/2)(1 + (2 + EPS/2)(1 + EPS/2))(\a®b\ + \c®d\)/nrm 

< (EPS)(\a®b\ © \c®d\)/nrm 

+ (EPS/2)(3 + 3EPS/2 + (EPS/2) 2 )(\a ®b\ + \c® d\)/nrm 

< (EPS)(1 + EPS/2)(\a®b\ © \c®d\)/nrm 

+ (EPS/2)(3 + 3EPS/2 + (EPS/2) 2 )(l + EPS/2)(\a ® b\ © |c ® d\)/nrm 

< (EPS/2)(1 + EPS/2)(2 + (3 + 3EPS/2 + (EPS/2) 2 )) (\a ® b\ © \c®d\)/nrm 
< (EPS/2)(1 + EPS/2)(5 + 3EPS/2 + (EPS/2) 2 )(\a ® b\ © |c ® d\)/nrm 

I 

We now use this lemma to get the error term for X/X. Of course, this error is less 
than the sum of the real and imaginary errors. Also, we let 

A = ((\x.re ® y.re\ © \x.im ® y.im\) © (\x.im ® y.re\ © |(— x.re) ® y.im\)) ® nrm 
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\(x .re<S>y .re@x .im®y .im)(Z>nrm— (x .rexy.re+x .imxy.im) / (y.rexy.re+y.imxy.im)\ + 
\(x .im<S>y .re(&(—x .re)<S>y .im)(Z)nrm—(x .imxy .re+(—x .re) xy.im)/(y.rexy.re+y.imxy.im)\ 

< (EPS '/2)(1 + EPS/2)(5 + 3EPS/2 + (EPS/2) 2 )((\x.re © y.re\ © \x.im © y.im\) 

+ (\x.im © y.re\ © \x.re © y.im\))/nrm 

< (EPS/2)(1 + EPS/2f(h + 3EPS/2 + (EPS/2) 2 )(((\x.re © y.re\ © \x.im © y.im\) 

© (\x.im © y.re\ © |(— x.re) © yim|)) nrm) 

= (EPS /2)(1 + EPS/2) 3 (5 + 3£PS/2 + (EPS/2) 2 ) (A) 
< (1 - (EPS/2)) 2 (5EPS/2)(1 + SEPS)(A) 
< (1 - (EPS/2))(5EPS/2)((1 + 3EPS) © A) 
< (5EPS/2) © ((1 + 3£PS) © A) 

Here we used the fact that 

(1 + EPS/2) 3 (5 + 3EPS/2 + (EPS/2) 2 ) < (1 - EPS/2) 2 (5)(1 + 3EPS) 

I 

Proposition 7.13 (A / A): 

x/y = (re,im;e) where 

re = (x.re © y.re © x.im © y.im) nrm where nrm = y.re © y.re © y.im © y.im 
im = (x.im © y.re x.re © y.im) nrm 

e = (1 + 4£PS) © (((5EPS/2) © A © (1 + 103EPS) © B) nrm) where 
A = (\x.re © y.re| © \x.im © yim|) © (\x.im © y.re| © |x.re © y.im\) 

B = x.e © (|y.re| © \y.im\) © (|x.re| © \x.im\) © y.e 

We are assuming (the computer will complain if this is not true) y.e < 100EPS © \y\, or, 
more accurately, 

y.e 2 < ((10000EPS) © EPS) © nrm. 

Proof: 

The error is bounded by \x/y — x y\ + (§)-e where 

, « | X | ~~ | •X' . & I I 

y \y\- y-e \y\ 
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\x\ 


\y\ 


+ 1 


y\ 


x.e 




x\ 


\y\ 


+ 


x\y.e 




\y\(\y\ 




y-e) 









\y\x.e + \x\y.e 

\y\(\y\ - y-e) 

< (1 + 101EPS) 



\y\x.e + 


x\y.e 




y\ 


2 



We have used the fact that our assertion y.e 2 < (lOOOOPPS EPS) nrm implies that 

' < (1 + 101EPS) ' 



\y\(\y\ - y-e) " |y| 2 ' 

which we derive in the next seven lines. 

y.e 2 < (10000PPS © EPS) nrm < (1 + (EPS/2))(10000EPS PPS) x nrm 

< (1 + (PPS/2)) 3 (10000PPS EPS) x (y.re 2 + yim 2 ) 

< (1 + (EPS/2)) 4 (10000EPS x PPS) x (y.re 2 + y.im 2 ) 

= (1 + (PPS/2)) 4 (100PPS) 2 x |y| 2 

This implies that y.e < A x |y|, where A = (1 + (PPS/2)) 2 (100PPS). Now, noting 
that j^j < 1 + 101PPS), we see that 

1 < i m i 1 ,, = r-r^- — it < (l + 101PPS) 1 



2 



|y|(|y|-y-e) |y|(|y| - A\y\) \y\ 2 {i-A)- ' \ y 

Resuming with our main proof, and setting 

A = (\x.re y.re\ © \x.im © y.im\) © (\x.im © y.re| © \x.re © y.im|) 
P = x.e © (|y.re| © \y.im\) © (|x.re| © |x.im|) © y.e 

we have that 

x 

\x/y - x y\ + (-).e 



< |:r/y - x y| + (1 + 101PPS) 



y 

'y|x.e + |x|y.e 



I?/ 



2 



< (EPS/2) (1 + EPS/2)(5 + 3EPS/2 + (EPS/2) 2 )((\x.re © y.re| © \x.im © y.zm|) 

+ (|x.im © y.re\ © |x.re © y.im\))/nrm 
+ (1 + 101PP5')(x.e x (|y.re| + \y.im\) + (\x.re\ + \x.im\) x y.e)/(|y| 2 ) 
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< (5EPS/2) (1 + EPS '/2) 2 (1 + (SEPS/10) + (EPS 2 / 20)) ((\x. re © y.re\ © \x.im © y.im| 

© (\x.im © j/.re| © |x.re © y.im\))/nrm 
+ (1 + 101£PS)(1 + EPS/2)(x.e x (|y.re| © |y.im|) 

+ (|x.re| © \x.im\) x y.e)/(( - ^gj^ ) 2 nrm) 



< (5EPS/2)(1 + EPS/2) 3 ((\x.re © y.re\ © |x.im © y.im\) 

© (|x.im © y.re| © \x.re © y.im\))/nrm 
+ (1 + 101PPS)(1 + PPS/2) 3 (1 + EPS/2) 2 (x.e © (|y.re| © |y.im|) 

© (|x.re| © \x.im\) ®y.e)jnrm 

= (1 + PPS/2) 3 ((5PPS/2)(A) + (1 + 101£PS)(1 + EPS/2) 2 (B))/nrm 
(1 + EPS/2) 3 ((5EPS/2)A + (1 + 103PPS)B)/nrm 
(1 + EPS/2) 6 (((5EPS/2) © A © (1 + 103EPS) © P) nrm) 
(1 + 4PPS) © (((5EPS/2) © A © (1 + 103EPS) © B) nrm) 
Here we have used the fact that 

(1 + EPS/2) 2 (1 + (3EPS/10) + (EPS 2 /20)) < (1 + (EPS/2) 3 ) 



Proposition 7.14 (V^l 7 ): 

Let s Q = y/ (\x.re\ © h a (x.re, x.im)) © 0.5 and <i = (x.im s)© 0.5, then = 
(re, im; e) where 

re = s if a;. re > 0.0 and re = d otherwise, 

im = d if x.re > 0.0 and im = s a otherwise, 

e = EPS © ((1 + AEPS) © (1.25 © s © 1.75 ®\d Q \)) 

Proof: 

This will be a little nasty. Let's begin by analyzing e s , which is the difference between 
the true calculation of s and the machine calculation of s, that is e s = \s — s \. First, we 
bound s. 

s = a/ (jaTrej + h(x.re, x.im)) * 0.5 
< (1 + EPS) 1/2 ^/(\x.re\ + h Q (x.re, x.im)) * 0.5 
< (1 + EPS) 1/2 (1 + EPS/2) 1/2 v / (\x.re\ © h a (x.re, x.im)) * 0.5 
< (1 + PPS) 1/2 (1 + PPS/2) 1/2 (1 + EPS/2) {J(\x.re\ © h Q (x.re, x.im)) * 0.5 

= (1 + EPS) 1/2 (1 + EPS/2) 3 ' 2 s 
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By a power series expansion, we see that 

(1 + EPS) 1/2 (1 + EPS/2) 3/2 

= (1 + \EPS - \ePS 2 + ...) + (1 + -EPS/2 + \{EPS/2) 2 + ...) 
2 8 2 8 

= (1 + ^EPS - ±EPS 2 + ...) + (1 + - A EPS + ^(EPS) 2 + ...) 
= (1 + ^EPS + ^EPS 2 + ...) 

So that, 

s<(l + ^EPS + ^EPS 2 + ...)s 

Similarly, 

s > (1 - ^EPS)s 
Thus, we can bound the s error, 

6 S — | S S Q | 

<((l + ^EPS+^EPS 2 + ...)-l)s 

= {\eps + ^eps 2 + ...) So 

Next, we analyze e^, which is the absolute value of the difference between the true 
calculation of d and the machine calculation of d. That is, ea = \d — d Q \. 

ed = \x.im/(2s) — x.im (2s )| 

< \x.im (2s Q ) — x.im/ (2s a ) \ + \x.im/ (2s a ) — x.im/ (2s) | 

< (EPS/2)\x.im/(2 So )\ + 

2 ss a 

T* 1 777 I 

< (EPS/2)\x.im/(2s )\ + {{h/A)EPS + (11/32)EPS 2 + ...)s a \ 

2 ss 

T 1 777 1 

< (EPS/2) \x.im/(2s ) | + | — — -{{5/4)EPS + (11/32)EPS 2 + ...) | 

2 s 

rp f TV! j 

< (EPS/2)\x.im/(2s )\ + I 2 So{1 - {5 /4 )E PS) {m)EPS+{UmEPS2 + - )l 
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< (EPS/2)(1 + EPS/2) _ { ljl ]Eps) \x.im (2»„)| 

= (EPS/2)(1 + EPS/2) (1 _ (5 7 / 4 2 )Epg) W 

Finally, we can bound the overall error e = e s + e^. 

e s + e d 

" + lk EPS2 + - )S ° + ( EPS ^ + EPS ^ \l-(I/4)EPS) ld ° l 

< + ^EP* + ...)(^ ) + ^(1 + ^) (1 _ (5/ 4 )gPg) l^l 

< SPS(1 + ^/2) (l _ (5/ ^ pg) (*,.) + EPS(1 + ^/2) (1 _ 1 ^| 

<m ( i + m/2 ) (1 _ (5/ ^ (^|la 

< EPS{1 + EPS/2f {i _ { ^ )Eps) { \ ® .„ © |I s 

< EPS(1- {EPS/2)){l + AEPS){^®s ®\ 7 -®d \) 

< EPS <g> ((1 + 4£PS) ®{^®s ®\-® d \)) 



Finally, we develop a couple of formulas for the absolute value of an XComplex. 

Formula 7.0 (absUB(X)): 

If x is an XComplex, then we get an upper bound on the absolute value of x as follows. 

\x\ = h(x.re, x.im) < (1 + EPS)h (x.re, x.im) 
< (1 - EPS/2)(1 + 2EPS)h Q (x.re,x.im) 
< (1 + 2£P5') <g> h (x.re, x.im) 

Thus, we define 

absUB(x) = (1 + 2EPS) <g> h a (x.re, x.im) 
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Formula 7.1 (absLB(X)): 

If x is an XComplex, then we get a lower bound on the absolute value of x as follows. 

\x\ = h(x.re, x.im) > (1 — EPS)h (x.re, x.im) 
> (1 + EPS/2)(1- 2EPS)h D (x.re,x.im) 
> (1 — 2EPS) © h (x.re, x.im) 

Thus, we define 

absLB(x) = (1 - 2EPS) © h Q (x.re, x.im) 

I 

Finally, in several places in the programs verify and fudging we perform a standard 
operation on a pair of doubles and must take into account round-off error. This is easy if 
we use Lemma 7.0. 

For example, in inequalityHolds we want to show that whxwh > absUB(along), where 
wh = absLB(whirle). By Lemma 7.0, we know that (1 — EPS) © (wh © wh) < wh x wh 
and we simply test that (1 — EPS) © (wh © wh) > abslIB (along). 

Similar situations occur in the functions horizon and larger-angle in fudging. 

A slightly more complicated version of this occurs in the computer calculation of 
pos[i] and size[i], that is, the center and size of a sub-box. Prior to multiplication by 
scale[i] = 2^ 5- ^/ 6 , the calculations of pos and size are exact. However, multiplication 
by scale introduces round-off error. For the center of the box we will have the computer 
use pos[i] © scaZe[i] with the realization that this is not necessarily pos[z] x scale[i}. Thus, 
we have to choose appropriate sizes to ensure that the machine sub-box contains the true 
sub-box. 

Notationally, this is annoying, because we typically use a computer command like 
pos[z] = pos[z] © scale[i], while in an exposition, we need to avoid that. We will denote the 
true center of the box by p[i] and the machine center of the box by po[i], and the true and 
machine sizes will be denoted s[i] and sq[i\. We will let pos[i] and size[i] be the position 
and size (true and machine are the same) before multiplication by scale[i}. 

Let p[i] = pos[i] x scale[i], po[i] = pos[i] <S> scale[i], and s[i] = size[i] x scale[i]. We 
must select sq[i] so that po[i] + so[i] > p[i] + s[i]. (Here, taking + on the left-hand side is 
correct, because the need for machine calculation there is incorporated at other points in 
the programs.) So, we must find sq[i] such that sq[i] > (p[i] — po[i]) + s[i]. 

(p[i] -po[i]) + s[i]. < (EPS/2)\p [i]\ + size[i] x scale[i] 

< (EPS/2)\p [i}\ + (l + EPS/2)(size\i] ® scale[i]) 

< (1 + EPS/2)((EPS/2)\p Q [i]\ + (size[i] ® scale[i})) 

< (1 + EPS/2) 2 ((EPS/2)\p Q [i] \ © (size[i] ® scale[i])) 

< (1 + 2EPS) ® ((EPS/2)\p [i}\ © (size[i] ® scale[i})) 

Thus we take 

s [i] = (1 + 2EPS) © ((EPS/2)\p [i]\ © (size[i] © scale[i])) 
This also works to give p [i] — s [i] < p[i] — s[i]. I 
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Chapter 8: Aff Approx's with Round-Off Error 

In Chapter 6, we saw how to do calculations with AffApprox's. In this chapter, we 
incorporate round-off error into these calculations. 

Convention 8.1: Recall that an AffApprox x is a five-tuple (x. f; x.fo, x.fi, x.f 2 ; x.err) 
consisting of four complex numbers (x.f,x.fo,x.fi,x.f 2 ) and one real number x.err. In 
Chapter 6, the real number was denoted x.e, but it seems preferable to use x.err in this 
Chapter. Also, we will suppress mention of the r functions used in Chapter 6 (see Definition 
6.4). As such, statements in this section will have to be translated (implicitly) to account 
for this suppression. 

Remark 8.2: One approach to round-off error for AffApprox's would be to replace the 
four complex numbers by four AComplex numbers complete with their round-off errors, 
and similarly for the one real number. We will not do this because it would necessitate 
keeping track of five separate round-off-error terms when we do AffApprox calculations. 

Instead, we will replace the four complex numbers by four XComplex numbers ("exact" 
complex numbers) and push all the round-off error into the .err term. In particular, in 
doing an AffApprox calculation, our subsidiary calculations will generally be on XComplex 
numbers and produce an AComplex number whose .e term will be plucked off and forced 
into the .err term of the final AffApprox. 

Convention 8.3: In what follows, we will use Basic Properties 7.0 and Lemmas 7.0 and 
7.1. Also, the Propositions in Chapter 6 will be utilized; as such, the numbering of the 
Propositions is the same in both Chapters (for example, Propositon 6.7 corresponds to 
Proposition 8.7). 

-X: 

Proposition 8.1: If x = (x.f; x.f , x.fi, x.f 2 ; x.err), then 

-x = (-x.f; -x.fo, -x.fi, -x.f 2 ; x.err) 



X + Y: 

We analyze the addition of the AffApproxs x = {x.f '; x.J 'q, x.J \, x.J ' 2 ; x.err) and y = 
(y-f',y-fo,y-fi,y-f2] y.err). To get the first term in x + y we add the XComplex numbers 
x.f and y.f; which produces the AComplex number r_f = x.f + y.f, and then we pluck 
off the XComplex part, r_f.z. The round-off error part r_f.e will be foisted into the overall 
error term r jerror for x + y. Similarly for the next three terms in x + y. 

Abstractly, the overall error term r jerror comes from adding the round-off error con- 
tributions r_/.e, r_/o.e, r_/i.e, r_f 2 .e and the AffApprox error contributions x.err, y.err. 
Of course, we have to produce a machine version. 

Proposition 8.2: If x = (x.f ;x.f , x.fi, x.f 2 ; x.err) and y = (y.f;y.f ,y.fi,y.f 2 ;y.err), 
then 

x + y = (rj.z; r.f .z, rji.z, r_f 2 .z; r -error) 
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where 

r-error = (1 + 3EPS) © ((x.e © y.e) © ((r_/.e © r_/ .e) © (r_/i.e © r_/ 2 .e))) 

Proof: The error is given by 

(x.e + y.e) + ((r_/.e + r_/ .e) + (r_/i.e + rj 2 .e)) 

< (1 + EPS/2)(x.e © y.e) + (1 + EPS/2)((rJ.e © r_/ .e) + (r_/i.e © r_/ 2 .e)) 
< (1 + EPS/2) 2 (x.e © y.e) + (1 + EPS/2) 2 ((rJ.e © r_/ -e) © (r^fi.e © r_/ 2 .e)) 

< (1 + £PS/2) 3 ((x.e © y.e) © ((r_/.e © r_/ .e) © (r_/i.e © r_/ 2 .e))) 

< (1 + 3£PS) © ((x.e © y.e) © ((r_/.e © r_/ .e) © (r_/i.e © r_/ 2 .e))) 
To get the last line we used Lemma 1 in Section 5. I 

X - Y: 

Proposition 8.3: If x = (x.f;x.f ,x.f 1 ,x.f 2 ; x.err) and y = (y.f; y./ , y./i, y./ 2 ; y.err), 
then 

x-y = (r.f.z; r.f .z, r.fa.z, r_f 2 -z; Terror) 

where 

r_/ = x.f -y.f 
r_f k = x.fk - y./fc 

r-error = (1 + SEPS) © ((x.e © y.e) © ((r_/.e © r_/ .e) © (r_/i.e © r_/ 2 .e))) 

I 

X + D 

Here, we add the AffApprox x = (x. f; x.fo,x.fi,x.f 2 ; x.err) to the double y. The 
only terms that change are the first and the last. 

Proposition 8.4: If x = (x.f; x./ , x./i, x./ 2 ; x.err) and y is a double, then 

x + y = (r-f.z; r_fo.z, r_f\.z, r_f 2 .z; r -error) 

where 

r_/ = x.f + y 
r_f k =x.f k 
r^error = (1 + EPS) © (x.err © r_/.e) 

Proof: The error is given by 

x.err + r_f.e 
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< (1 + EPS) <g) (x.err © rj.e) 
by Lemma in Section 5. I 
X - D 

Proposition 8.5: If x = (x.f; x.fo, x.fi, x.f 2 ; x.err) and y is a double, then 

x — y = (r-f.z; r-fo.z, r-f±.z, r_f 2 .z; r -error) 

where 

r_f = x.f -y 

r-fk =x.fk 
r_f.z = x.f -y 
r^error = (1 + EPS) © (x.err © r_/.e) 

I 

X * Y 

We multiply the AffApproxs x, y while pushing all error into the .err term. 
We will use the functions (see Formulas 7.0 and 7.1, at the end of Chapter 7) absUB = 
(1 + 2EPS) © hypot Q (x.re, x.im) and absLB(x) = (1 — 2EPS) © hypot (x.re,x.im). 
When x is an AffApprox, we define dist(x) to be 

(1 + 2EPS) © (absUB(x.fo) © (absUB(x.fi) © absUB(x.f 2 ))). 

This is the machine representation of the sum of the absolute values of the linear terms in 
the AffApprox x (the proof is straightforward) . 

Proposition 8.6: If x = (x.f; x.f , x./i, x.f 2; x.err) and y = (y.f;y.f ,y.f 1 ,y.f 2 ;y.err), 
then 

x * y = (r-f.z; r_f .z, r_f\.z, r_f 2 .z; r -error) 

where 

r_f = x.f x y.f 
r_f k = x.f x y.f k + x.fk x y.f 
Then, x x y = (r_f.z; r_f$.z, r_f\.z, r_f 2 .z; r -error), where 

r -error = (1 + 3EPS) © (A © (B © C)) 

with 

A = (dist(x) © x.e) © (dist(y) © y.e) 
B = absUB(x.f) © y.e © abslIB(y.f) © x.e 
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C = (rj.e © r_f .e) © (rj 1 .e © rj 2 .e) 

Proof: We add the non-round-off error term for x x y to the various round-off error terms 
that accumulated. 

((dist(x) + x.e) x (dist(y) + y.e)) + ((absUB(x.f) x y.e 

+ absUB(y.f) x x.e) + (rj.e + r_f .e) + (rji.e + rj 2 .e)) 

< (1 + EPS/2) 2 [(dist(x) © x.e) x (dist(y) © y.e)} + (1 + EPS/2)(absUB(x.f) © y.e 

+ absUB(y.f) © x.e) + (1 + EPS/2)((rj.e © r_/ .e) + (r_/i.e © r_/ 2 .e)) 

< (1 + EPS/2) 3 [(dist(x) © x.e) © (dist(y) © y.e)} + (1 + EPS/2) 2 {(absUB(x.f) © y.e 

© absUB(y.f) © x.e) + ((r_/.e © r_/ .e) © (r_/i.e © r_/ 2 .e))} 

< (1 + EPS/2) 3 A + (1 + EPS/2) 2 (B + C) 

< (1 + EPS/2) 3 A + (1 + EPS/2) 3 (B + C) 

< (1 + EPS/2) 4 (A®(B + C)) 

< (1 + 3EPS)®(A®(B + C)) 

I 

X * D: 

Proposition 8.7: If x = (x.f;x.fo,x.fi,x.f 2 ',x.err) and y is a double, then 

x x y = (rj.z; rj .z, rj±.z, rj 2 .z; r_error) 

where 

r_/ = x.f xy 
rj k = x.f k x y 

r-error = (1 + 3EPS) © ((x.e © |y|) © ((rj.e © r_/ .e) © (r_/i.e © r_/ 2 .e))) 

Proof: 

(x.e x |y|) + ((r_/.e + r_/ .e) + (r_/i.e + r_/ 2 .e)) 
< (1 + EPS/2)(x.e © + (1 + EPS/2)((r_f.e © r_/ .e) + (r_/i.e © r_/ 2 .e)) 
< (1 + EPS/2) 2 (x.e © |y|) + (1 + EPS/2) 2 ((r_f.e © r_/ .e) © (r_/i.e © r_/ 2 .e)) 

< (1 + £PS/2) 3 ((x.e © |y|) © ((r_/.e © r_/ .e) © (r_/i.e © r_/ 2 .e))) 

< (1 + 3EPS) © ((x.e © |y|) © ((rj.e © r_/ .e) © (rjx.e © r_/ 2 .e))) 
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I 

X/Y: 

For convenience, let ax = absUB(x.f), ay = absLB(y.f). 

Proposition 8.8: If x = (x.f;x.f ,x.f 1 ,x.f 2 ;x.err) and y = (y.f;y.f ,y.f 1 ,y.f 2 ;y.err), 
then 

x/y = (r-f.z; r_fo.z, r^fi.z, r_f 2 .z; r^error) 

where 

r_f = x.f/y.f 
rj k = (x.f k x y.f - x.f x y.f k )/(y.f x y.f) 
r.error = (1 + 3EPS) (((1 + 3EPS) A (1 - 3EPS) B) © C) 

with 

A = (ax © (dist(x) © x.e)) B 

B = (ax ay © dist(x) ay) © ((dist(y) ax) (ay ay)) 

C = (r_/.e © r_/ .e) © (r_/i.e © r_/ 2 .e) 

B = ay (1 + EPS) (dist(y) © y.e) 

Of course, we do have to be concerned about division by zero, so the program will complain 
if ay is not greater than (1 + EPS) (dist(y) © y.e), that is, if D is not greater than zero. 

Proof: 

As usual, we add the round-off errors to the old AffApprox error, taking into account 
round-off error. Let's work on it bit by bit. 

(ax + dist(x) + x.e) / (ay — (dist(y) + y.e)) 

< (ax + (1 + EPS/2)(dist(x) © x.e)) /(ay - (1 + EPS) (dist(y) © y.e)) 

< (1 + EPS/2) 2 (ax © (dist(x) © x.e)) /(ay - (1 + EPS) (dist(y) © y.e)) 

< (1 + EPS/2f(ax © (eftst(x) © x.e))/(- 1 )(ay © (1 + EPS) (dist(y) © y.e)) 

1 © EPS 1 2, 

< (1 + BPS/2)4(ax © (dist(x) © x.e)) (ay 0(1 + BPS) (dist(y) © y.e)) 

= (1 + BPS/2) 4 A 
< (1 + 3EPS) A 
The next term, being subtracted, requires opposite inequalities. 

(ax j ay + dist(x) / ay) + dist(y) x ax /(ay x ay) 
> (1 - EPS/2) (ax 0ay + dist(x) ay) + (1 - EPS/2) (dist(y) ax)/( - _ ^ ps ^ ) (ay ay) 
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> ((1 - EPS/2) 2 (ax ay © dist(x) ay) + (1 - EPS/2) 3 (dist(y) ax) (ay ay) 
> (1 - EPS/2) 3 [(ax ay © dzs£(x) ay) + (dist(y) ax) (ay ay)) 
> ((1 - EPS/2) 4 ((ax ay © a'zst(x) ay) © ((dist(y) ax) (ay ay))) 

> (1 + EPS /2)(l + 3EPS)(B) 
> (1 - S 

Finally, we do the round-off terms. 

((r_/.e + r_f .e) + (r_/i.e + r_/ 2 .e)) 

< (1 + EPS/2)((rJ.e © r_/ .e) + (r_/i.e © r_/ 2 .e)) 

< (1 + EPS/2) 2 ((rJ.e © r_/ .e) © (r_/i.e © r_/ 2 .e)) 

= (1 + EPS/2) 2 C 
Now, we put these three pieces together. 

(ax + dist(x) + x.e) / (ay — (dist(y) + y.e)) 
— ((ax I ay + dist(x) / ay) +dist(y) x ax/(ay x ay)) + ((r_/.e + r_/o.e) + (r_/i.e + r_f 2 .e)) 

< (1 + 3£PS) A - (1 - 3£PS) B + (1 + EPS/2) 2 C 
< (1 + EPS/2)((\ + 3EPS) A © (1 - 3£PS) S) + (1 + EPS/2) 2 C 

< (1 + £PS/2) 2 (((1 + A © (1 - 3PPS) P) + C) 

< (1 + PPS/2) 3 (((1 + 3PPS) A © (1 - 3PPS) B) © C) 

< (1 + 3PPS) (((1 + 3PPS) A © (1 - 3PPS) B) © C) 

I 

D/X: 

We are dividing a double x by an AffApprox y. For convenience, let ax = |x|, ay = 
absLB(y.f). 

Proposition 8.9: If x is a double, and y = (y.f;y.fo,y.fi,y.f2;y.err), then 

x/y = (r_/.z; r_/o.z, r.fr.z, r.f 2 .z; r_error) 

where 

r_/ = x/y.f 
rj k = -(x x y.f k )/(y.f x y./) 
r _ error = (l + 3EPS) (((1 + 2EPS) (ax B) (1 - 3EPS) B) © C) 
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B = ax ay © (dist(y) ax (ay ay)) 

C = (rj.e © r_f .e) © (r_/i.e © r_/ 2 .e) 

D = ay © (1 + .EPS) (dist(y) © y.e) 

Again, we have to be concerned about division by zero, so the program will complain if ay 
is not greater than (1 + EPS) (dist(y) © y.e), that is if D is not greater than zero. 

Proof: 

Let's start with the pieces. 

ax/(ay- (dist(y) +y.e)) 
< ax/ {ay - (1 + EPS) (dist(y) © y.e)) 

< ax /( 1 + Ep S /2 ^ ay e (1 + EPS) ® ( dist ^ ® y- e ^ 

< (1 + EPS/2) 2 ax .D 
< (1 + 2£PS) (ax D) 
The next term, being subtracted, requires opposite inequalities. 

ax j ay + dist(y) x ax/ (ay x ay) 

> (1 - EPS/2)ax ay + (1 - EPS /2)dist(y) ax/( — ^ ps ^ )( a y ® ay) 

> (1 - EPS/2)ax ay + (1 - EPS/2) 3 dist(y) ax (ay ay) 

> (1 - EPS/2) 3 (ax 0ay + dist(y) ax (ay ay)) 

> (1 - PP5'/2) 4 [ax ay © (dist(y) ax (ay ay))] 

> (1 + EPS/2) (1 - 3EPS)B 
> (1 - 3£PS) £ 
Finally, we do the round-off terms as before. 

((r_/.e + r_/ .e) + (r_/i.e + r_/ 2 .e)) 

< (1 + EPS/2) 2 C 
Putting the pieces together, we have: 

< ax/ (ay — (dist(y) + y.e)) 

— (ax / ay + dist(y) x ax /(ay x ay)) + ((r_/.e + r_/o.e) + (r_fi.e + r_f 2 .e)) 

< (1 + 2EPS) (ax 0D)-(1- 3EPS) S + (1 + EPS/2) 2 C 

72 



HOMOTOPY HYPERBOLIC 3-MANIFOLDS ARE HYPERBOLIC 

< (1 + EPS/2)((\ + 2EPS) (ax D) (1 - 3£PS) S) + (1 + EPS/2) 2 C 

< (1 + £PS/2) 2 (((1 + (ax D) © (1 - B) + C) 

< (1 + PPS/2) 3 (((1 + 2PPS) (ax D) 9 (1 - 3EPS) B) © C) 

< (1 + 3PPS) (((1 + 2PPS) (ax D) © (1 - 3PPS) B) © C) 

I 

X/D: 

We are dividing an AffApprox x by a double y and the computer will object if y = 0. 
Proposition 8.10: If x = (x./; x./o, x./i, X./2; x.err) and y is a double, then 

x/y = (r_f.z; r_fo.z, r_fi.z, r_f 2 -£; r .error) 

where 

r_/ = x.f/y 
r_f k = x.f k /y 

r -error = (1 + 3PPS) ((x.e |y|) © [(r_/.e © r_/ .e) © (r_/i.e © r_/ 2 .e)]) 

Proof: 

This is easy. 

s-e/|y| + (( r -/-e + r_/ .e) + (r_/i.e + r_/ 2 .e)) 
< (1 + EPS/2)x.e |y| + (1 + EPS/2) 2 [(r_f.e © r_/ .e) © (r_/i.e © r_/ 2 .e)] 

< (1 + £PS/2) 2 (x.e |y| + [(r_/.e © r_/ .e) © (r_/i.e © r_/ 2 .e)]) 

< (1 + EPS/2f(x.e |y| © [(r_/.e © r_/ .e) © (r_/i.e © r_/ 2 .e)]) 

< (1 + 3EPS) (x.e \y\ © [(r_/.e © r_/ .e) © (r_f 1 .e © r_/ 2 .e)]) 

I 

Here, x is an AffApprox and we let ax — absUB(x.f). There are two cases to consider 
depending on whether or not D = ax (1 + EPS) (dist(x) © x.e) is or is not greater 
than zero. 

Proposition 8.11a: If x = (x./; x./o, x./i, x./ 2 ; x.err) and D is not greater than zero, 
then we use the crude over-estimate y/x = (0; 0, 0, 0; (1 + 2EPS) {J (ax © (xdist © x.e)) ) 

Proof: 

\J ax + xdist + x.e 
< \Jax + (1 + EPS/2)(xdist © x.e) 
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< + EPS/2) 2 (ax © (xdist © x.e)) 
= (1 + EPS/ 2) a/ (ax © (xdist © x.e)) 

< (1 + EPS/2) 2 ^/ (ax® (xdist® x.e)) 

< (1 + 2£PS) © v/(ax© (xc&st ©x.e)) 

I 

Proposition 8.11b: If x = (x.f;x.fo,x.fi,x.f2',x.err) and D is greater than zero, then, 

= (r_/.z; r_/o.z, r./i.z, r_/ 2 -^; r -error) 

where 

r_/ = a/x7/ 

t = r_f + r_f 

rj k = AComplex(x.f k .re,x.fk.im,0)/t 

(Simply put, r_fk = x.fk/(2y / x.f). The reason we have to fuss to define r_fk is because 
\Jx.f is an AComplex.) 

r -error = (1 + 3EPS) © ( 

{(1 + EPS) © Zfax 
© (1 - 3£PS) © [dist(x) (2 x $ax) ®</D}} 
© ((r_/.e © rjo-e) © (r_/i.e © r_/ 2 .e)) 

) 

Proof: 

Let's work on the pieces. 

y'ax < (1 + EPS) © ^ox 

Next, 

dist(x) / (2yfax) + \J ax — (dist(x) + x.e) 
> (1 - EPS/2) 2 dist(x) (2^ax) + ^ax - (1 + EPS) © (dist(x) © x.e) 

> (1 - EPS/2fdist(x) (2^/az) 

+ (1 - EPS/2) 1 / 2 v / ax © (1 + EPS) © (dist(x) © x.e) 

> (1 - EPS/2) 2 dist(x) + (1 - EPS/2) 3/2 v r D 

> (1 - £PS/2) 2 [cfet(x) (2^ax) + v^D ] 

> (1 - £PS/2) 3 [cfet(x) (2^ax) © v^D ] 
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> (1 + EPS/2){1 - 3EPS)[dist(x) (2^ax) © </D ] 
> (1 - 3EPS) [dist(x) © <[D ] 

Adding in the usual term, we get as our error bound 



\fax — (dist(x) I (2y/ax) + \/ ax — (dist(x) + x.e) ) + ((r_f.e + r_f .e) + (r_/i.e + r_/ 2 .e)) 



< (1 + £PS) <fax 

- (1 - 3£PS) [dist(x) (2^ax) © v^D ] 

+ (1 + EPS/2) 2 ((rJ.e © r_/ .e) © (r_/i.e © r_/ 2 .e)) 



< (1 + EPS/2){(1 + EPS) ^ 

© (1 - 3£PS) [dist(x) © v^D ]} 

+ (1 + EPS/2) 2 ((r_f.e © r_/ .e) © (r_/i.e © r_/ 2 .e)) 

< (1 + EPS/2) 3 ({(1 + EPS) ^Hx 

(1 - 3EPS) [dist(x) (2 -Vox) ® <fD }} 

© ((r_/.e © r_/ .e) © (r_/i.e © r_/ 2 .e))) 



< (1 + 3£PS) ({(1 + EPS) ^ax 

(1 - 3£PS) [dist(x) 

© <fD }} © ((r_/.e © r_f .e) © (r_/i.e © r_/ 2 .e))) 
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